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A tneoretical study of the time-dependent and steady 
hypersonic viscous interaction with surface mass transfer has 
been described in this work. The problems analysed here are 
pertinent to the two-dimensional hypersonic flow over a 
slender wedge-wing commonly employed as lifting surface in 
hypersonic vehicles. For the purpose of analysis, the 
transformed flow governing equations have been obtained for 
the general case of a time-dependent flow over a slender wedge 
involving the surface mass transfer . 

First the steady problem with homogeneous (same gas) 
and heterogeneous (foreign gas) injections is analysed. The 
surface injection presents a possible means of alleviating the 
aerodynamic heating problem which assumes formidable dimensions 
at hypersonic speeds. The use of injection is also made for 



XX 


controlling the boundary layer and inducing the control forces 
on aerodynamic vehicles. Most of the earlier works connected 
with the study of surface mass transfer are without the study 
of the corresponding viscous- inviscid interaction phenomenon . 
Some of the strong-interaction flows involving surface mass 
transfer have been treated under the restrictive category of 
similar solutions. The present study treats the non-similar 
problem in an exact way for the entire viscous-inviscid 
interaction regime (from 'strong ' to 'weak' including the 'transi- 
tion') by using the complete tangent-wedge approximation for 
the calculation of interaction-induced pressure in the layer 
contained between the shock wave and the boundary layer, A 
series solution has been developed to initiate the solutions 
at the beginning of the strong-interaction regime. The 
laminar boundary layer equations for a binary mixture of 
perfect gases have been numerically solved for surface 
injections of helium and argon with variable fluid properties 
and various values of the wall injection parameter. The 
solutions have also been obtained with the in jection/suction 
of air at the surface. The computed results include the 
conventional boundary layer parameters of displacement 
thickness, wall shear function, enthalpy gradient at the wall 
and viscous-inviscid interaction induced pressure. Tho 
comparative effectiveness of the various coolants for the 
high speed problem from the point of view of controlling the 
induced pressures is pointed out. The analysis of this steady 
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problem has been carried out by the two commonly employed 
numerical schemes due to Clutter & Smith and Marvin & 

Sheaffcr . The Clutter-Smith technique employing the shooting 
method is highly accurate but fails to converge for certain 
problems with surface injection where the shin-friction and 
heat transfer become small. The finite-difference method 
(using the Crank-Nicolson scheme)due to Marvin & Sheaffer, 
although somewhat less accurate, has the distinct advantage of 
employing the dependent variables rather than their slopes to 
obtain solutions. Therefore , this method could provide results 
for the cases where the Clutter-Smith technique did not converge. 
The conparative advantages and disadvantages of the two methods 
are discussed in detail for the problems involving viscous- 
inviscid interactions. 

The results from the steady problem have been further 
utilised to obtain the lift and drag coefficients, the lift 
to drag ratio, the position of the aerodynamic centre, and the 
heat transfer rate per unit width of the wedge-wing. 

Next the time-dependent problem with surface injection/ 
suction of air has been analysed. The problems of maneuvera- 
bility and control of hypersonic vehicles require the study of 
unsteady flows such as the one over a slender wedge after an 
impulsive change in the angle of attack of the wedge-wing. 

In the present work, the time history of such a time-dependent 
flow, till the final steady state is reached, is analysed and 



xx ii 


tile time -dependent solutions of the flow governing equations 
with air in jection/suction are used for determining the 
temporal variations in the pressure, wall shear and heat transfer 

rate numerically. Time-dependent problem with foreign gas injection 
has not been considered here due to the large computational time 
required for such problems. 

The parameters employed for defining both the steady 
as well as unsteady problems are the free stream Mach number 
M . Reynolds number R based on the free stream conditions 

00 0 T ft 

°o,li 

and the length of the wedge 1'*, the semi-wedge angle 0^ and 
the surface injection/suction parameter, a. For the unsteady 
problem, in addition, the initial and final angles of attack 
a\ and oc f are also defined. For the problems considered 
here and R >> 1 and the angles involved 0-, , a - , 

<< 1. However, the relative magnitudes of these various 
parameters are such that the combination fC/(R^ )V2 

00 jL 

remains finite and the product M^j 6 [ £0(1), where e is the 
maximum deflection of the wedge surface. Finally, two typos of 
thermal wedge conditions on the wedge surface are used for the 
steady case: insulated wedge surface and the wedge surface 
maintained at a constant temperature. For the unsteady problem, 
however, only the latter thermal condition is employed. A 
comparison with the available results, wherever possible, is 
included. 



CHAPTER I 


STATEMENT OP THE PROELM AND ASSOCIATED LITERATURE 

1.1 Introduction 

The two problems analysed here are the problems of 
steady and unsteady hypersonic flow over a slender wedge-wing 
commonly employed as lifting surface in hypersonic vehicles. 

The surface heating of these wedge-wings assumes formidable 
proportions at hypersonic speed and the means are to be 
provided to control this heating. Cooling by the homogeneous 
or heterogeneous surface injection in the boundary layer 
provides an effective way of doing it. The behaviour of 
wedge-wings in steady hypersonic flow situations without the 
surface-injection control has been studied in detail during 
the past three decades. Not much attention, however, seems 
to have been given to the problem involving the surface injection/ 
suction. Further, the problems of maneuverability and control 
of hypersonic vehicles require the study of unsteady flow such 
as the one over a slender wedge after an impulsive change in 
the angle of attack of the wedge-wing. We are particularly 
interested in the transient behaviours of the pressure 
distribution, shear stress and heat transfer rate along the 
sides of the wedge, which will furnish informations about the 
lift, drag and other aerodynamic coefficients. Thus, the 
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significance of the problems analysed here lies in determining 
the control characteristics and devising the effective ways 
of controlling the surface heating of hypersonic vehicles 
using wedge-shaped aerofoils. 

Before discussing the associated literature with the 
problems outlined here, some more insight is provided into 
the unsteady flow phenomenon and the surface mass transfer. 
Taking the unsteady problem first, Pig .1.1 » schematic ally 
represents such a flow over a slender wedge-wing. 



Pig, 1.1 Plow over a slender wedge-wing before and after the 
impulsive change in angle of attack. 


The governing parameters are of the following order 
of magnitude; 


and 


>> 1 * 0 ( 10 ) 
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where p* , u* , y* denote density, velocity and viscosity, 
respectively, in the free stream and L* is the length of 
the wedge . 

Since for a slender wedge with a sharp leading edge 
the shock is attached to the leading edge, the flow on one 
side of the wedge is independent of the flow on the other side 
and, therefore, we can consider them separately. The problem,, 
therefore, reduces to determining the unsteady hypersonic 
viscous flow over an inclined surface following an impulsive 
change in the angle of inclination. The disturbances due to 
the presence of the wedge are confined within a thin curved 
shock be.yond which the free stream conditions prevail. Similar 
model of an inclined surface (with the shock attached at the 
leading edge) replacing a wedge-wing may also be used for the 
steady problem. Thus, the flow-model shown in Figure 1.2 
for the hypersonic flow over an inclined surface (or flat 
plate) may be used for analysing the problems under 
consideration • 

In any unsteady viscous flow problem there are two 
separate time scales involved % 

(i) ~ characteristic time given by L*/ u * , where 
L* is a representative body dimension and u* is 
the free stream velocity . 

(ii) t* v - characteristic time for the diffusion of heat 
and vorticity in the boundary layer. 
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EDGE 

1. Free-nolacular flow region 

2. Continuum- flow region with slip-boundary condition 

3. Strong- inter act ion region 

4. Transition region 

5. Weak- interaction region 


Fig 


1.2 Flow Regions Over an Inclined Flat Plate 
in Hypersonic Flow (Schematic) 
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In hypersonic flow these two time scales are of the same 
order of magnitude i.e* 0(10“^). Let T be the time in which 
a change in the flow conditions or the geometry of the solid 
body occurs. A nonsteady problem is that in which T is of 
the same order of magnitude as or m nv or less. By 

impulsive change we mean a change occuring in a time T 

<< T iny , In the other extreme case i*e. T >> X j Lllv .» T v 

the flow may be considered as quasi- steady . 

I he unsteady as well as steady hypersonic viscous flow 
over a slender body is characterized by the presence of a high 
temperature, low-density boundary layer adjacent to the body. 

This boundary layer is much thicker than the corresponding 
boundary layers in subsonic or moderately supersonic flows* 

As a result of the curvature of this thick boundary layer, 
a sufficiently large outward streamline deflection is induced 
to bring about an important change in the " effective ’ ' 
geometric shape of the body, with the result that pressure 
variations are propagated into the external inviscid layer. 

This external pressure field in turn feeds back into the boundary 
layer and affects its rate of growth. Thus the boundary layer 
and surrounding inviscid flow are mutually interdependent. 

This phenomenon, referred to in the literature as hypersonic 
viscous interaction or shock-boundary layer interaction, 
therefore, requires the solutions for a steady or unsteady 
boundary layer subjected to an * ’ external (streamwise) pressure 
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gradient’ ' which, is not known a priori, but which depends on 
the rate of growth of the boundary layer itself through the 
hypersonic inviscid flow relations* 

Coming to the problem of surface injection, this presents 
a possible means of alleviating the formidable aerodynamic 
heating problem as pointed out earlier. In addition, the use 
of injection is also made for controlling the boundary layer 
and inducing the control forces on aerodynamic vehicle. The 
injection of the same or different gas through the boundary 
layer materially influences the rate of heat transfer to the 
surface, skin-friction and the viscous-interaction induced 
pressure field. Under the hypersonic flow conditions, the 
pressures induced by the thick boundary layer are greatly 
increased with surface mass transfer. For cooling purposes, 
either the injection of the same gas or that of a foreign gas 
(e.g, helium) may be used. If the injected gas is at the 
surface temperature, it would act as an ' ‘insulant’’ by 
altering the velocity profiles so that the heat conducted to 
the surface is reduced and, accordingly, the surface temperature 
would be smaller. In addition, the diffusion of a foreign 
gas would also reduce the viscous dissipation. 

A review of literature on the interaction problem in 
two-dimensional steady/unsteady hypersonic viscous flows with 
or without surface mass transfer is presented in the following 
section. 
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1.2 Review of Literature 

The flow regions, shown in Pig. 1.2, on an inclined 
flat plate are according to the model introduced hy Lees and 
Probstein and contained in Q,2j , In the basic flow model 
postulated by them, two main flow regions can be distinguished; 
I, a ''leading-edge'' region? II, a region of no- slip . 

In the leading edge region I, two sub-regions exist* 

The relatively small domain ’’l'' around the leading edge is 
of the order of 20 mean free paths in breadth, and the concepts 
of Maxwell and Smoluceowski cannot be applied. This is also 
known as the free-molecular flow region. Downstream of 
zone 1 the lateral extent of the viscous layer is sufficiently 
large compared with the mean free path so that the Maxwell- 
Smoluchowski concept of a small region of slip and temperature 
jump at the surface with a main body of continuum flow is 
probably valid. Here, the shock emerges from the viscous 
region, and a distinct, larger inviscid layer exists. Zone 2 
may be defined as the slip-flow region. 

In region II, slip and temperature jump at the surface 
are, generally, negligible, and all concepts of continuum flow 
apply throughout. The domain of no-slip is divided into three 
sub-regions; (3) ''strong'' interaction region? (4) transition 
region? (5) ''weak'' interaction region. As we proceed 
downstream along the wedge the induced effects decay, and 
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eventually if the Reynolds number, Re becomes sufficiently 

°°,x 

high these effects must be weak. In the weak- interaction 
region (zone 5) the effects produced by the self-induced 
pressure gradient are essentially perturbations on the classical 
high Reynolds number flow. In the transition region (zone 4), 
as the name suggests, the interaction is neither ' ' strong' 1 
nor ’ ’weak* ! . 

In part of region I and completely in region II, the 
flow field between the shock and the body surface is divided 
into two distinct regions, the external inviscid layer and the 
viscous boundary layer. In the external layer, the flow is 
governed by the inviscid Euler equations, whereas in the other 
region, the flow is governed by the boundary layer equations [[321 
as long as (6*/ x *) << 1. Thus, the hypersonic viscous 

interaction problem, mentioned earlier and as applicable to 
region II, has the following three aspects: 

(i) The flow in the external layer can be considered 

as the inviscid hypersonic flow over an * ' effective body 1 * given 
by the boundary layer displacement thickness added to the 
original body . 

(ii) The distributions of pressure, tangential velocity 
and total enthalpy at the edge of the 1 ' effective body 11 
determine the flow in the boundary layer . 

(iii) The displacement thickness used to obtain the 
external inviscid flow should be consistent with the 
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displacement thickness obtained from the solution of the 
boundary layer equations subjected to the external pressure 
gradient . 

The steady-flow interaction problem appears to have been 
studied in detail. Various approaches used to analyse the 
interaction problem for the steady flow are reviewed in detail 
by Hays and Probstein £2^] , Dorrance J~4J , Stewartson J”b3 
and Moore []6^j . Since our interest in the steady-flow problem 
is mainly under the conditions of surface mass transfer, only 
the basic characteristics of these approaches are outlined 
here briefly . 

The flow in the external layer, which is equivalent to 

the inviscid hypersonic flow over a sharp-edged body of specified 

shape can be obtained by the method of characteristics, the 

shock expansion theory and the tangent -w edge approximation. 

The last method is the easiest of the three, gives results of 

good accuracy and has been used extensively in the hypersonic 

interaction problems. In the tangent wedge approximation the 

pressure at any point on a slender body of thickness y*(x*) 

6 

is approximated by the pressure across an oblique shock that 
produces the local flow deflection dy*/dx*. Por free stream 
Mach number M >> 1 and for a slender body with dy*/dx* << 1, 
with the finite product of K(x*) = M^dy */dx*, the tangent-wedge 

G 

approximation may be written as DO: 
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For small and large values of K, equation (1.1) can be 
expanded as: 
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For the hypersonic interaction problem K in the equation (1.1), 
(1*4) or (1,5) has to be replaced by: 



d6*(x*) 
dx^ - 1 


( 1 . 6 ) 


where y*(x*) defines the body surface and 6*(x*) is the 
boundary layer displacement thickness. Since the displacement 
thickness 6* in (1.6) is not known a priori , the use of 
tangent wedge approximation still retains the basic iterative 
nature of the interaction problem. However, the advantage 
obtained is that the inviscid flow in the external layer need 
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not be calculated explicitly. A disadvantage of employing 
the tangent wedge approximation is that the vorticity in the 
external layer due to the shock curvature, which may affect 
the structure of the boundary layer, can not be accounted for. 
However, for a sharp-edged slender body, the effects of the 
vorticity interaction can be neglected £23 beyond a region 
very close to the leading edge. 


The extent of the interaction between the boundary 
layer and the inviscid external flow is characterized by the 
hypersonic interaction parameter, x , defined as: 


with 
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( 1 . 8 ) 


Here C is the coefficient of the linear viscosity-temperature 
relation. 


For the strong-interaction flows, x >> 1 and eq. (1.5) 
is the appropriate expansion of the tangent-wedge approximation. 
For the weak interaction flows, eq. (1.4) is to be used 
corresponding to X < 0(1). For the general case, where a 
continuous solution is sought from strong to weak interaction 
region through the transition region, no simplifications to 
eq. (1.1) can be made. 
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The foregoing discussion gives some insight into the 
interaction problem in steady hypersonic flows without the 
surface mass transfer. With homogeneous (same) or heterogeneous 
(different) surface mass transfer, the interaction problem 
becomes little more complex. Most of the earlier works 
connected with the analysis of surface mass transfer are without 
the study of the corresponding interaction problem. We will 
briefly describe in the next few paragraphs some of these works 
pertinent to the surface mass transfer problem. 

The basic features of surface mass transfer effects are 
found in the early studies of the low-speed boundary layer. 
Schlichting do and Emmons and Leigh ^8] investigated the 
incompressible similarity boundary layer on a flat plate with 
air injection at the surface. Their results indicate marked 
reduction in wall shear. Compressibility effects were first 
included in studies of air injection by Lees [[9]] , Eckert and 
Livingwood Q.(7j , Low [If] , and others. Lees investigated the 
stability of the compressible laminar boundary layer with air 
injection, Eckert and Livingwood carried out simplified 
calculation to obtain the effects of air injection and 
comparisons of convection, transpiration, and film coolimgj 
and Low completed the exact solution for compressible flow 
with v fe « x 1//2 and constant Prandtl number. All investigators 
confirmed the benefits of air injection upon the reduction of 
heat transfer and equillibrium wall temperature. 
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In a recent work, Inger and Swean jjl£] have presented 
a comprehensive study of similarity solutions for vectored 
injection and suction of air in non-adiabatic laminar boundary 
layers with zero axial pressure gradient* However, in a 
hypersonic flow with surface mass transfer, it is necessary 
to consider the self induced pressure gradients resulting 
from the displacement of the external flow field due to the 
presence of a relatively thick boundary layer. 

The effects of vectored air injection in the * 'weak- to- 
moderate* * interaction region have been investigated by 
Inger and Swean jjL3] . They treated the interaction effect 
as a small, non-similar perturbation caused by the self-induced 
pressure gradient on a basic self-similar isobar ic flow 
solution with heat transfer and vectored air injection. 
Considering both upstream and downstream vectored injection, 
their results show that the physical properties of the flow 
situation are significantly influenced by both the velocity 
and angle of injection. Zien [14] investigated the unsteady 
boundary layer on a flat plate with surface mass transfer. 

Results on skin-friction, heat transfer rate and weak interaction 
pressures were obtained under conditions of surface mass transfer 
that correspond to similarity solution in steady flow. 

The foregoing discussion was on cases of homogeneous 
injection in which the external flow and the injected fluid 
are of the same medium. When the injected fluid is different, 



14 


a binary mixture boundary layer flow results with the diffusion 
phenomenon having an important influence. 

In a detailed investigation, Baron jjL5] derived the 
complete equations for laminar flow of a binary mixture on 
the basis of the equations of change resulting from Enskog's 
solution of the Boltzmann equation. It was shown that although 
the mass and momentum conservation equation remained unchanged, 
additional terms appeared In the energy equation to account 
for the interaction between energy and mass fluxes. The 
descriptive system was augmented by a component mass-continuity 
relation, The equations were solved to yield similar solutions 
for a flat plate with zero axial pressure gradient, both 
Prandtl and Schmidt numbers being taken as unity. Helium and 
Carbon-di-oxide injections were considered. Greater reductions 
in skin-friction and heat transfer were reported for the 
injection of gases lighter than air as compared to heavier gas 
injection. Wuest Q.5] obtained similar solutions with variable 
fluid properties whose dependence on temperature and concentration 
was represented by semi-empirical formulae. Integral forms 
of the boundary layer equations were also derived. A more 
general analysis of the problem (not restricted to similar flows) 
was made by Moran and Scott [if] • The equations were solved 
in terms of Crocco variables using a finite-difference technique. 
The method requires starting profiles and is not apparently 
capable of handling large changes taking place in the streamwise 
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direction. Culick [18] presented an approximate analysis of 
the compressible turbulent boundary layer on a porous flat plate 
with distributed surface mass transfer, Li and Gross jj-9] 
obtained similar solutions for the hypersonic laminar boundary 
layer on a flat plate considering strong-interaction and 
normal surface injection of a foreign gas. A set of analytical 
expressions were derived relating the induced pressure at the 
wall, wall heat transfer, and skin-friction in terms of 
parameters describing the strong interaction flow. However, 
the mathematical restriction on body shape and injection 
velocity necessary for the existence of similar solutions 
disallow the treatment of problems of practical interest, 

Jaffe, Lind and Smith jjjoj investigated the hypersonic 
laminar boundary layer on a cone with surface mass transfer » 
and described a general method for solving such a problem, 
numerical results were obtained for injection of helium, argon 
and air, but the hypersonic interaction phenomenon was not 
considered. Jain and Li [21] investigated the binary boundary 
layer flow characterised by finite external stream Mach number 
and finite pressure gradient. A Gor tier’s type series 
solution was developed. They carried out detailed computations 
of the thermodynamic and transport properties of hydrogen-air 
and helium-air mixtures and presented a number of graphs to 
illustrate their behaviour. 
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Blottner [22} solved the Boundary layer equations for a 
multicomponent flow with finite chemical reactions. He 
evolved a scheme for solving the governing equations for a 
flow problem involving a large number of chemical species and 
obtained results in the case of a cone and a hyperboloid 
subjected to a high Mach number flow. Marvin and Sheaffer [23] 
suggested a numerical method for solving non-similar and similar 
boundary layer equations, including binary gas injection for 
non-reacting ga„ses. They got the starting profiles by solving 
the similarity equations obtained by putting x=0 in the 
non-similar case. This approach is not mathematically sound 
since it assumes the validity of the boundary layer equations 
at the leading edge. They did not analyse the associated 
interaction problem also. 

Gupta and co-worhers [24,25] very recently presented 
results for the case of homogeneous and heterogeneous vectored 
injection cooling at hypersonic speeds. They also obtained 
the first order correction to the displacement induced pressures 
for the problem of strong-interaction flows which may not be 
treated under the category of similar solutions. However, 
the complete interaction problem involving the recomputation 
of the flow field from the modified induced pressure field 
was not considered. 

How, coming to the unsteady flows, the interaction 
problem, in comparison with the steady state, is further 
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complicated because the lateral velocity of the unknown 

t 

displacement surface also affects the external flow. Lighthill 
[20 presented an extension of the tangent wedge approximation 
to obtain the time-dependent pressure on a slender moving body 
of shape y*(x*,t*) placed in an inviscid hypersonic stream. 

This extension was based on the 1 'piston analogy’ ' of Hays [20 
who observed the equivalence between the steady two-dimensional 
flow over an arbitrary body and the unsteady one- dimensional 
flow in front of a moving piston. The pressure in front of 
a piston moving with a velocity w* in a stationary compressible 
medium is given by 
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where a* is the velocity of sound and the subscript °° refers 
to conditions in the undisturbed gas. Equation (1.9) is analogue 
to the corresponding tangent wedge approximation (1.1) if one 
identifies w* with the normal velocity component induced 
because of the slope of the body and neglects quantities of 
0(M" ) . Eor a moving body an additional term 9y*/dt* is 
introduced in the normal velocity component. Thus the 
instantaneous pressure p*(x*,t*) on the body is given by 
eq. (1.9) with 
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Milos [28] has discussed the expressions (1*9) and (1.10) 
and its application to the unsteady pressure in inviscid 
hypersonic flow. For the small and large values of K, the 
expansions of expression (1.9) are similar to equations (1.4) 
and (1.5) with the meaning of K given by eq. (1.10). 'The 
first two terms of the small K expansion, equation (1.4), 
lead to the acoustic approximation formula, while the piston 
theory of Lighthill [ 26 ] corresponds to the first four terms. 

For the unsteady interaction problem y* in eq. (1,10) 
should be replaced by (y* + &*), where A*(x*,t*) is the 
displacement thickness for the unsteady boundary layer, 
introduced by Moore and Ostrach [29] . 

Examples of unsteady interaction problems include the 
works of Reshotko, Rodkiewicz, Gupta & co-workers [30,31,32,33, 
They considered the unsteady weak- inter action flow on a semi- 
infinite flat plate at zero- incidence after the free stream 
Mach number was impulsively increased/decrcased by a small 
amount (of 0(1$ )). The pressure on the plate ’’was assumed 
constant for calculating the unsteady boundary layer and a 
small unsteady perturbation was calculated from the boundary 
layer solutions. This approach is justified because of the 
small change involved in the problem. Gupta and Rodkiewicz 
[34,35,36] later considered the unsteady strong interaction 
flow for the stepwise-accolerated flat plate in the direction 
of its original motion. 
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The problem of impulsive motion of a flat plate involving 
the development of the flow from initial unsteady state described 
by Rayleigh to the ultimate steady state described by Blasius 
has been solved by Hall j~37] . His results were limited to the 
case of an incompressible flow. Hall employed the natural 
coordinates, namely, the three independent variables in space 
and time. Stewartson [38,39J and Dennis and Walker [40, 4l] 
obtained solutions in terms of the two transformed independent 
variables. Their results are identical to those of Hall [37] . 
Ban and Kuerti [45] have also made important contributions in 
the area of unsteady flows with particular application to shock 
tube flow. Gupta [43] extended their work to the case of a 
non-linear problem. These works are now overtaken by those 
of Dennis and Stewartson. 

Based on the work of Marvin and Sheaffer Nath 

[44,46] and his co-workers treated the unsteady compressible 
boundary layer flow. Their results, however, do not include 
the viscous interaction effects. 

The unsteady problem involving the binary-diffusion of 
perfect gases was analysed by Gupta and co-workers [46,47] 
with application to expansion tube flows. This analysis 
without the binary-diffusion, was also extended to the case of 
unsteady turbulent boundary layer flow [48,49] . 

The time-dependent interaction problems involving the 
lateral motion of the body are solely confined to cases of 
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harmonic oscillations. King C50J and Or 1 ik-Ru ck emann j__5 13 
considered the interaction on ah oscillating slender wedge. 

These are examples of quasi-steady analysis, that is the flow 
at any instant was assumed to be the steady f low corresponding 
to the conditions prevailing at that instant. Demetri and 
Gupta |_ 52 U k av e elao analysed the compressible oscillating 
boundary layers but their results are limited to supersonic 
flow problems without incorporating the viscous interaction 
effects . 

file interaction problems involving impulsive lateral 
motion of tne body are very few in the literature £533 . The 
method of solution is based on ’ quasilinearization ’ which 
essentially is an extension of the IMewt on-Rap hson method of 
solving tne algebraic and transcendental equations. The results 
of ref. [354] do not have a mono tonic ally increasing variation 
with time. However, the wavy behavior of the induced pressure 
and tne gradient of the enthalpy at the wall is nor explained 
and there is no physical basis to justify this wave-like 
variation in time. 

1.3 Description of the Present Work 

The present work is concerned with the study of the 
steady and unsteady high speed viscous- inter act ion f lows which 
may not be treated under the category of similar solutions. 

This analysis has been undertaken with the objective of evaluating 
the viscous interact ion- induced pressures, displacement thickness, 
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shear-stress and wall heat-transfer functions etc. over a 
slender wedge-wing or equivalently an inclined flat plate 
moving at hypersonic speeds with surface injections of air, 
argon and nelium. For the steady flow case, variable fluid 
properties are used and the solutions are obtained for various 
values of the wall injection fluxes. fne viscous interaction 
pressure is calculated by employing the tangent-wedge approxima- 
tion and the variables employed are similar to those of refs. 

[53] and For the unsteady problem involving the small 

impulsive lateral motion of a slender wedge-wing, an extension 
of the tangent-wedge approximation is used. For this case, 
however, the constant fluid properties are used with the surface 
injection of air only. 

In Chapter II of the present work, the details of the 
derivation of transformed governing equations are given. Chapter 
III deals with the steady flow over an inclined flat plate with 
air injection/suction at the surface. A four- term series 
solution needed to supply the starting profiles (and other 
details) for the strong-interaction region is developed here. 

This solution is used to specify the Initial conditions in 
space to carry out the numerical integration from strong to the 
weak interaction regime. An elaborate (and more accurate) 
numerical integration scheme (sometimes known as "difference- 
differential" or Clutter-Smith technique) employing the 
shooting method is described in the latter parts of this chapter; 
and the results are discussed toward its end. Chapter IV deals 
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with the details of the numerical integration for the above 
method. 

Chapter V is concerned with the heterogeneous injection 
(employing argon and helium) under the steady flow conditions. 

In the beginning portions of this chapter the shooting method 
for argon injection is presented; and in the letter part a 
f inite-dif ference scheme after envisaging a new transformation 
is described, Tne main results are also discussed at the end 
of this chapter. 

Cnapter VI indicates the utilisation of results, of the 
previous chapters, on the two sides of the wedge yielding the 
characteristics of the two-dimensional wing for steady flows 
with and without surface mass transfer. Similar results for 
the unsteady problem are not given due to the prohibitive 
computer time required for such calculations. 

last Cnapter deals with the unsteady problem with air 
injection/suction on an inclined flat plate following an 
impulsive change in the single of inclination. The required 
initial conditions in time are derived from the steady state 
solutions of Chapter III and the solution to the unsteady problem 
is sought by a finite-difference scheme similar to the one used 
in Chapter V. Here the finite-difference scheme is preferred 
over the "difference-differential” technique employing shooting 
method for the reasons explained in this chapter. 
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THE FLOW GOVERNING EQUATIONS 


2.1 Time-Dependent Two-Dimensional Compressible Boundary 
Layer Equations for a Binary Gas-Mixture of Perfect 
Gases 


The unsteady laminar boundary layer equations for a 
binary gas-mixture with no chemical reactions are considered 
for the analysis of two-dimensional high speed flow over an 
inclined surface. Tne governing equations in the Cartesian 
Co-ordinate system with origin fixed to the leading edge of the 
fla.t plate are Hl5] : 
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Continuity of Species Equation : 
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Equation of State : 

p* = R* p* T* 

And the Viscosity-temperature-species Concentration Relation ; 
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Thus, the non-dimensional form of the equations (2,1) through 
(2,6) e,re : 
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6p 


0 


ay 


(2.9) 


( 2 . 10 ) 


OH dH OH 0 1 OP 

at ax + ay H p at 


[ H -u 2 1 





1 + c^(s - 1) 
1 + c 1 (r - 1) 


"f 


Re 


~,L* P L 


fy % §> + % < ~ } + W { * (Le " 1} 


(• 


■) (H-u ) > 


-f-~f 


Q + 

r/ & _ 

lu p — 

9c- 9c. Qc,. 

5 “ + U — + V ay“ 

and 


Re. 




1 8_ ( le ffij 

p ay \Pr ay ' 


(2.11) 


( 2 . 12 ) 


p = 5 pf 


u - y(e^jl) 


(2.13) 

(2.14) 


m A 

m = 2 , c - 

m i 

e = m*/mj 

P = 


where 
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lie 


°°,L* 


u p L 

OO oo 


* 


y* C* 


Pr 


k 


Le 


n # -n* 
p °12 ° d 

k* 


The details of simplifications employed in obtaining equations 
(2.9) and (2.11) are given in Appendix A. 

2.2 Transformations of the Governing Equations : 

The governing equations are now transformed to bring 
them to a form suitable for numerical solution. First we 
introduce a stream function, f , defined by 


PU 


djj 

ay 


(2.15) 


and the Dorodnitsyn-Howarth transformation 


x = x 


J 

y = J p dy 

^b 


(2.16a) 

(2.16b) 


t = t (2.16c) 

Now for any function c(x,y,t) we obtain the following relations 


dx 


it 

ay 


= Is + is 

dx dy 


dx 

P is. 

ay 


(2.17a) 


(2.17b) 
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Lz = ££ SiL + &£ 
dt dt 3y at 

Then Equations (2.15) and (2.8) reduces to 

dip 


(2.17c) 


u = 


V = 


ay 


1 (M + §2 M + 2Z) 

0 aS 8x 0 y 6t 


(2.18) 

(2.19) 


With the help of relations (2.16) through (2.19), equations 
(2.9), (2.11) and (2.12) can be written as i 

1 + c 1 (s-l)' 


<LjL. + M 9Jl_ _ M = _ | 2 [h - (£i) 

3yat 3y 3x8y dx dy 2 2 p dx LI ay 


1 + c^(r-l) i 


1 £_ / x c n £l) 

e ~,L ay ay 


(2.20) 


aH + aji an __ a_£_ 


at ay ax dx ay 

+ ~ — ~ i - — ( 


a=E 


pH - (M)21 


” 1 + c 1 (s-l)“ 

L ay -J 

i 

__ 1 + c 1 (r~l) _ 


Re 


“,1* 


a , xc PaHx.„a , * c p " 2 


Pr — 

ay ay 


2) + 2 £- (Pr-1) ^ s-| ) 

” Qy rJ - 


34 o 

ay ay 5 


+ (Le-lJ (-- fel ) Xi , 

ay r c r x i ^ i ay ay 


3c. 


( 2 . 21 ) 


!fi+ 

ay ax 


94 9c i 


ax ay ^,l* ay v Pr 


ft x o p ac. 

n Le - 1 


) 


ay 


( 2 . 22 ) 


where 


yp 
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For homogeneous injection 



v = C I 

(2.23) 


p = pi 

(2.24) 

and 

yp = G p 

(2.25) 


where C is the Chapman-Rubesin linear viscosity constant and 
is given by 


C 


jil/2 

1 b 


1 + S*/£* “1 

' CO 

- \ + s*aL 


(2.26) 


The values of S* in the above relation is taken as 110°K. 


With heterogeneous injection 


up = X G p (2.27) 

where x = — — and has a value of unity for air injection or no 
o p 

injection case. At this stage we introduce the normalised 
pressure function 

p(x,t) = — -r- - (2.28) 

p 0 X 

and change x, y, t to x, n , t defined by 



x = x 

(2.29a) 


^ ^ f Re co,L* 1 l/2 

~ - TT — 172- ) 

2(p 0 x x) ' 

(2.29b) 

and 

t = t 

(2.29c) 
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Boting that x = X L * 


wn ere 


Z = M 3 {_£ ) 1 / 2 

X L* u oo ^Re J 


(2.30) 


°°,L* 


is the interaction parameter evaluated at the edge of the 
wedge, we can write 


A y x 


1 

"4 


where 


1 , Re <» t Ij* ^ 1 /2 

2 


A = * 


(2.31) 


(2.32) 


C P 


o XL* 


is a constant, for any function ?(x,y,t) we have 


further 


5c _ 8 c n 3 c 
ax ax 4x on 

1 


~ = A x 4 

ay 

a c _ S£ 
at at 

we define 


Si 

3n 


f (x, 


n ,t ) = A 



t|> (x ,y , t ) 


(2.33a) 


(2.33b) 

(2.33c) 


or, ^ = f (2.34) 

A x 

Now using expressions (2.33) and (2,34) in (2.20), (2.21) and 
(2.22) we get 
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& s!i_ + 4S (“ kn . « 2^) . f 4 . p a_ (X 

StQn n 3xQn 3x 6n p 6n^ n 3 n p 


+ p & % 

p 3x 



r ,31,21 


~ 1 + c 1 (s-l)~ 

1) 

h - fc) 


i 


1 on 

L J 


L 1 + c 1 (r-l) _ 


= 0 (2.35) 


45 <£ + 45 (®£ as . |S H) . t as 

d^E an a5 8,1 as 


3 n 


x 

- -=• 

p 3t 




r i + c 1 (s-i)“ 

I . 

__ 1 + c 1 (r-l)__ 

i - p 

i 


3n W an^ + 2 5n f ibc ( Pr ~ 1 ) 


3f off 
d“ Qn 2 


- P 


r a 


f - r-1 


9n 


b (le-l) { 


Pr 


c- (r-1 )+l 


} 


{H 


/Qfx 2 ^ C 1 

k an ' 1 ar~ 


= 0 


(2.36) 


= d°i _ 

4x FT + 4x (fs 


af 8o i 


3x 


.p S c * Oo* -v \ Qc« 

M —i) - f i - p (-— Le '—■) 
Q = 3n J 9n ^ 3n v Pr 3 n 


0 


(2.37) 

The co-ordinates x and t have not been distorted in 
the transformations. This facilitates the retrieval of 
physically significant results from the solutions of the 
transformed equations. The distortion of the co-ordinate 
normal to the surface provides a uniform computational domain 
as the boundary layer thickness varies from a very small value 
(close to the leading edge) to a fairly large value at the 
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weak- interact ion region. Thus we can write equations (2.35) 
through (2.37) in the form : 


4x 


p 

S f 
0 tOn 


+ 4x ( 


Of 8 2 f 
On OxO^ 


dff Of \ 
On 2 9X 


a 2 f - 

• — ^ - p 

e n 


0__ 

On 


( X 


2 

Of 
3 n^ 



H 


- (^) 2 ' 

K br\ J 


1 + C -j^( S— 1 ) 

1 + c 1 (r-l) 


0 (2.38) 


4 X SH 4 /Of. OH 

^x at + (, Qn Qx 


OH Of 
d n Qx ' On 


) - f H - 2 If 


H 


(M)2 

''On' 


“ 1 + C^(s-l) ~ 

- 0 

X 

_ 1 + c ^ ( r- 1 ) _ 

“ p On | 

1 

[_Pr 


|f + 8(Pr-U §£ 2-| 


+ (Le-1 ) { 


r-1 


c 1 (r-l)+l 


} { H - (— ) 2 > 

on 9 n 


= 0 (2.39) 


Qc, 


4x 


8T + fa 1 


3c- ~..p 3c - 3c . p, ^^4 

~ i. - i.) - f ~ - p (^- Le rr^) = ° 

On Ox Ox: On ' On x 9n v Hr On 


(2.40) 


where p = p(x,t) is given by the expression (2.51) given in next 
section. 


In the next section we discribe the equations required for 
interaction with the flow in the external layer. 


2.3 The Interaction Equations : 

For the unsteady flow over a flat plate inclined at 
a small angle 0^ to the main flow, the ’’effective body” shape 
is given by 
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y*U*,t*) = e b x* + &*(x*,t*) 


(2.41) 


Tile normal velocity component at the edge of this 
effective body is given by equation (1.10) 


w*(x*,t*) = u* I 9. + in 

5x • 


(2.42) 


Here a*(x ,t*) is the displacement thickness of the unsteady 


boundary layer, given by Moore and Ostrach 


At high 


speeds the displacement and boundary layer thickness are 
approximately same. With the present notation the differential 
equation for a* becomes 


P e U e A *'A (P*u* - p*u*) dy< 




(2.43) 


Using the continuity for the external flow 


^ + ^( p * u * ) = 0 


(2.44) 


equation (2.43) may be simplified to 


e at e ax 


p* a* / (1 

6 e * v 


—if) ay* 

P u 

e e 


jL a V P 

p e L ay 

y-u p. 


(2.45) 
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Applying the transformations of the last section on the two 


integrals appearing equation (2.45) it is seen that 


* 

/ 6 (1 
y? 


* * J e 

-f~) ay* * /. (i 

p e u e K 


•) dy' 


(2.46) 


(y — 1 ) 

viq M „ 


3/4 °° 1 + Ch(s-I) 

r 


p o L 1 + c 1 (r-1) 


Ii - (4-) 4n 
v 5 n 


(2.47) 

_2 

where terms of 0(11^) have been neglected from the integrand of 
the last integral. Writing 6* for both the integrals in 
equation (2.45) and using (2.44) once again, we get 


5 A* . * 9 A* _ Q6* . * d6* 

+ a e — - — + u e — 

6 ij 5 x 5 v ox 


(2.48) 


A further simplification is obtained if u* is replaced 
x 0 

by u*. This is consistent with Stewartson cd and Goldsworthy 
[56] that 


e * 

u 

00 


u _2 

0(M ) 


finally we get from (2.42) and (2.48), after non- 


d imen s i on ali s ing , 


£(x,t) 


/ fl , 66 , 66^ 
^b + ax + 


(2.49 ) 


<-* ^ 1 Y x T * 3/ 4 00 r i+c/s-Dir o .p o *" 

c / o y- 1 * li* x ' r 1 u / 0 i \2 ^, n 

6 (x,t ) = — = — - — -J — 7 — 7 H - fc) dn 

I fp^ M « p o L l+c 1 (r-l) J [ 0t) j 


(2.50) 
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where p = 


P 0 x 


Po x 


1 + YK 


Y ±i K + {(Iti K) 2 + 1) 



L * J 

For flows involving surface mass transfer, displacement 
tnickness A* or tne boundary layer thickness 6* snould be 
replaced by their effective values |~ 57~j s 


( 2 . 51 ) 


A * 

A eff 


x* pJ v* 


A_ ^ = A + J — r dx 


ef f 


6 * + J" 


x 


* * 
p e u e 

p iW 

*: < 


(2,52) 


dx* 


* 


(2.53) 


-.vhere a* and 6* are the values without surface mass transfer. 

file last terms of the above equations are due to the 
injected mass flow at the surface. How from the equation of 
state (2.6), for the quantities across the boundary layer we 
may wr ite ; 

p: 


b. 

* 


* 


1 


(2.54) 

'e “b 

From the hypersonic assumption of equation (1.3) we also have 

H* 


1 — ^2 ^ 1 i m P- L yi n S that I* % 0. More appropriate, 


7T 11 

2 oo 


it is implied that I*/T* « 1. Therefore, the contribution from 
the surface mass transfer effect may be neglected in equations 
(2.52) and (2.53) at hypersonic speeds. 


(xdd&d -fcY Coirr? ^>6- (l *>-£££ 
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2.4 Boundary and Initial Conditions 

The boundary conditions in n are given by the 
conditions on the wedge surface and at the outer edge of the 
boundary layer. On the surface of the wedge the streamwise 
component of velocity is zero and the wedge is either insulated 
or held at a constant temperature, Thus, 




P 

* 

ii 

o 

(2.55a) 

Y* = V*( X *) 

(2.55b ) 

°i = °id x n 

(2.55c) 

— = 0 (for insulated wedge surface) 

ay 

(2.55d ) 

I* = I* = constant (for constant temperature 

wedge surface) 

(2.55e) 


Bor a given problem, v*(x*) and c-j^x*) are related 
to each other. Eckert has suggested that the convective and 
molecular transport terms for air must be balanced at the surface 
with foreign gas injection. Thus, for a zero net flux of the 
free-stream gas, namely, air at the surface, Eckert proposed 
the following relation between v*(x*) and c^Cx*) for a 
physically realistic solution from the aerodynamic point of 


( P*V*) 


b 




b 


v iew j 
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In terms of the transformed variables we have 


at n = 0 


Sf 

QtT 


f + 4x || = - p, v, (x) 


2 V fie oo t L* Tf 


dx ^ b -,rrr” r-~ p* 




Sc 


sr - <7Te> “ * < f + 4 * i> 


(2.56a) 

(2,56b) 

(2.56c) 


SH 

8n 


0 


or 


Mxj = 4 




(2.56d) 


b VA ' ^b 
1 u *2 ~ 
2 


r c lb + ( 1 — c 11^) 


lb' 


21, 


(y-1) K 


(2.56e) 


where 1, is the constant wall temperature arid r = c* / c* • 

P 1 p 2 

At the edge of the boundary layer the tangential velocity and 
the total enthalpy should match with the corresponding Quantities 
in the external layer and the injected species concentration 
vanishes. Thus 


at y* = y* 


u* = u* 
e 


H* = H* 
e 


c^ = 0 


(2.57a) 

(2.57b) 

(2.57c) 


In terms of the transformed variables we have 
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at n “ 

§f = 1 (2.58a) 

H = 1 (2.58b) 

c 1 = 0 (2.58c) 

-2 

where "terms of order have been neglected. 

In addition to the above conditions, in the normal 
direction, the initial conditions in x and t are also required. 
These are discussed in chapters III, V and VII. 

In boundary layer calculations it is customary to 
satisfy the conditions (2.58) at a finite value n = n . The 
value of is chosen sufficiently large such that in the 

resulting computations Jf«( Tl ) - ij < 10 and |H(n )-l| <_ 10"“ 4 . 

Prom experimental executions a value of n = 5.5 is found to 
be sufficient to suffice this condition. 



CHAPTER III 


THE SI'S ARY FLOW OS AN INCLINED FLAP PLATE WITH HOMOGENEOUS 

INJECTION 


3.1 Introduction 

Equations (2.38) through (2.40) may be specialised for 
the steady hypersonic flow on a flat plate inclined at a small 
angle 0^ to the free stream and with air injection on the 
flat surface. With air injection the equation of state takes 
the form given by equation (2.24) and the viscosity law is 
given by equation (2.23). They are 


p = P T 


and 


= c T 


Thus the momentum equation (2.38) and energy equation (2.39) 
become 


P 


si 

8n^ 


* f rf - - & H } - * 


» 3f af a f 


8 n 


8 n* 


2x 

p dx 


tj / 8f\2 


(3.1) 


and 

L. of? 1 4. f 4y /9f 5H 8H af x 

Pr Qn 2 8n k 8n 8x “ 8n dx' 


L 8n dn 8n 3 


Pr “v — ^ Qn 2 
For this case p = p(x) is given by (2.51) with 


(3.2) 



and 6 = 6(x) given by (2.50). Equations (3.1) and (3.2) are 
the same as those given in ref, [58) for the no-injection case. 

3.2 Boundary and Initial Conditions 


Ihe boundary conditions for this problem may be obtained 
by taking c lb (x) = 0 in equations (2.56) and (2.58). Then 
equation (2.56c) can be dropped and (2.56e) becomes: 

H = - constant 

Eor obtaining the initial condition in x, a series 
expansion solution is given in the next section. 

3.3 Series Expansion Solution 

In this section a scheme for solving the steady strong 
interaction case by series expansion is described. Making use 
of proper expansions QQ for p,6,f and H in terms of the 
interaction parameter x and k b (= e b ) the partial differenti- 
al equations (3.1) and (3.2) are reduced to a sequence of 
ordinary differential equations. The unknown constants 
appearing in the expansion of p,6, etc, are then determined 
from the solution of these ordinary differential equations in 
accordance with the tangent-wedge equation. In this section 
solutions for terms up to third order are carried out for 
arbitrary but a constant value of Pr, 

Por the steady problem p = p(x) is given by the tangent- 
wedge approximation of equation (2.51) : 
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P = 


P r> X 


1+YK 


Y+ 1 Y + 1 \ ^ 1/2 

K + { ( T -ti K) + 1 } 


with K(x) = M (6,+ 

00 b dx 


and 

constant 


6 = 6(x) is given by (2.50). If we define a 


lc b 


k, = M 0, 
b oo b 


then 


d6 


K(x) = k. + M ™ 
v ' b «> dx 


(3.4) 


(3.5) 


Por large values of IC(x) the above expression for p gives 
(see equation (1.5;) 


P = 


PnX 


lil+ll k 2 + K“ 2 + 0(K“ 4 ) 


y-b 1 


(y+1) 1 


Por the strong interaction region (x >> 1) an order of 
magnitude analysis PQ shows that 


or 


M - oc 
00 x 


I, I ^ « 
°° dx 


- V2 

X 

(3.6) 

-1/2 

X 

(3.7) 


1/2 


Writing K(x) = k^ + constant x and substituting in the 
expansion (1.5) for p for large value of K, we have for , 
x » 1 


p(x) = P 0 x 


1 + !A + kiA. + 0(f 3/2 )' 

X X 


(3.8) 
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where p^(i=0,l ,2,3) are constants to be determined. 
Consistent with (3.8) we have for 6 


6(x) _ ^o x 
x ” M 


_!/2 


i + + + o(x' 3/2 )' 

•*** 7 \/ 


(3.9) 


where 6^(i=0,l,2,3) are constants. 

In view of the expansions (3.8) and (3.9), Lees and 
Probstein suggested the following expansions for f and H: 

(n , + alf % + LiLLLLf / 

o 

and 


i.tlk ...-3/2. 

~~T 7 ' g ~ + a 

X 


f(x,n) - f„(n) + — — t 7 "q " + — " — 3f__ — + 0(x ) (3.10) 


H (n)k H (n)+H„(n)k h ^-3/2 

H(x,n) = H o (n) + -—-^- 7 ^'- + 3__b_ + 0 (x ) (3.U) 

X * 

Using the expansion (3.10), we may also write 

2 


f , !_ at _ f (n) , f l (n)k b , V n) + V n > V , 

t + 4X ~ - f 0 (n) + -“Tyg" + 1 — “ + 


f l'(n)k b 

+ - i T7V + 2 


f (n) + f,(n) k 


•b 


+ . 


Hence from (2.56b) for the steady problem we can obtain: 

1/2 


- PV 


so that 


Re T * 1/2 

°°,1 x 7 


T CT 


o 


3172 ' 2 


f , , , „ , . f 2 U > + f 3 <">V , 

f 0 ( n) + 2 • - YJ2 + ^ ~~Z " — — — + . 


Re 


V 


■ p b v b 


°°i2_ 


r. 


Vc P 0 X L 


8f 1 (o)t b fgCoMjtOt,. 

f (0) + — + 3 — +■ 


T 72 ' 


X 


(3.12a) 
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The equation (3.12a) gives the required, form of the 
surface mass flux PgV^ so that for x ■* 00 similarity structure 
of the solution is obtained. This, inject ion/ suet ion is 
performed according to the relationship : 


p b v b Re 


1/2 


.x’ 


fC p - 1 / 2 
V k o X 


— oc + 


a 0 a l 

+ -- + 


X 


where a, a Q , a^,... are all constants of order unity. 
If we consider tne special case a n = 0 for n = 0,1,2, 


f 0 ( 0) = -2a 


f x ( 0) = fg( 0 ) = f 3 (o) 

iVe can also rewrite (3.1 2b) as 


0 


v b = 


r ~ 0 


- 1/2 

X 


V" Re * 

V oo,X^ 


a 


a + 


o 

Jj2 

X 


+ a i. 


-1 


j 


(3.12b) 


. . , then 
(3.13a) 
(3.13b) 


(3.12c ) 


a = (v- b ) Q with a n = 0 (3.12d) 

Thus, in seeking the series solution, v^ should be of tne form 
given by equations (3.12c) or (3.12d). However, a . form 

of v^, used with the numerical solution of non-similar equations 
(3.1) and (3.2) can be of the type given below i 


* _ 

L Pb . 


v b ^b^o + R 1 — ^ (3.12c) 

With Re # v 0(10 ) the normal velocity in equation (3.12-0) or 
(3.l2d) will be small relative to the tangential velocity as 
long as a v 0(1). Tne value of x at tne beginning of the strong 
interaction region is ^0(10 ). lor this value of x, it is noted 
that (6/x) « 1. Thus, the boundary layer approximations are 
valid in the hypersonic flow as long as this condition is met J/74]. 
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where is a constant of 0(10“ ). In terms of the dimensionless 

stream function eq* (3.12e) may be written as 


(f + 


af 


2K- 


4x ^i). = -2a ±-~ 

Sx b 


][We 


“jX 1 


*1‘ 


LfT 


(3 ,13a) 


or, after replacing the density by pressure with the help of 
the equation of state, 


(f + 4x 


5f 

5x 


) = -2a 


4K fp M p 
1 v ^o °° 

(r-D H b ~ 



(3.13b) 


Next, the x-derivatives of 6 can be obtained from (3.9). 
Then from (3.5) we can get K(x). Using this K(x) in (1*5) 
and then comparing with (3.8), we can obtain the following 
relations between the two sets of constants 6^(i=0,l,2,3) and 
P ^ ( i= 0»^-»2,3) • 


p 0 = Is < 1,+1) 6 c 2 

(3.14a) 

si = 1 (6 i + jy 1 

(3.14b) 

10 c .32 3Y+1 1 

p 2 3 2 9 y(Y+ 1) § b 0 2 

(3.14c) 

p 3 = 6 3 + M (5l + i-) 2 

( 3 *14d ) 


On inserting the expansions (3.8), (3.10) and (3,11) in the 
equations (3.1) and (3.2) and equating the terms of the same 
order on the two sides we get the following sets of coupled 
ordinary differential equations for f^(n), H^( n) (i=0,l,2,3): 
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! * 1 


+ Vo + = 0 


h" + Prf H* - 2(1-Pr ) 

0 0 O v ' 


" 11,2 i Ht 

<V - Vo 


= 0 


t i 


! I 


f l + f o f l - £ o £ l + 2 Vl + P 
= P 1 j~2 ^V f o ^ ~ f o f o J 

f| f ! t Iff If M III I s 

H x + Prf 0 H 1 _Erf 0 H 1 -S(l-Er)(f 0 f 1 +2f Q fj, +f Q fj.) 
= -a?V £ i + p i Erf o H o 

4” + £ o £ 2 "<*2 + 3f o f 2 + P<W2 f 


= P. 


2 P ( H „- f o 2 > + Vo 


] 


«• t I N , f Iff If ft Ml 1 

H 0 + Pr£JH 0 - 2Prf Ji 0 -2(1-Pr) (fj n +2f Q f 2 +f Q f ; 


'2 ‘ * ““0**2 ““cr2 

= “ 3M o f 2 + *S? rf o H !) 


o 2 


f 3 +f o f 3 “ 2f o f 3 + 3f o f 3 + ^ H 3" 2f o f 3^ 

■ p 3 [ + Vo] - Pi I (V £ o e ) 

+ ^1 2 ^ H l“ 2f O f l' ) ~ ^1 3 -2 ^1^1 + ^ + P 


t t 


H 3 o H 3 - 2Er£ o H S + 3PrH o f 3 

I Iff II If S II ! . 

- 2(1-Pr) (f 0 f 3 + 2f 0 f s + f 0 f 3 ) 

I I Ilf M It HI t 

= P 3 Prf 0 H 0+ p 1 j2(l-Pr)(Vl + 2£ o £ 1 +£ o V “ 

fff ft P. t I . 

+ 2(1-Pr) j^f 1 f 1 + (f 1 ) j+ Pr(f 1 H 1 -2f 1 H 1 ) 


(3.15a) 

(3.15b) 

(3.16a) 

(3.16b) 

(3.17a) 

J 

(3.17b) 

(3.18a) 

II — 

i] 


(3.18b) 
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Here primes denote differentiation with respect to n . In 
equations (3.16) through (3.18) the unknown constants 
p^(i = 1,2,3) occur in the inhomogeneous parts. To remove 
these unknowns we introduce a further change of variables 
following ref. [64] 


1 

_ "-<r» 

««. ii 

H, = — 

(3.19) 

Pi 

1 p l 

f 0 



2 

P2 ’ 

H 2 = pf 

(3.20) 

II 

CnI 

3 f 2 )/ Pl 2 , H 3 = (H s - p 3 h 2 )/ Pi 2 

(3.21) 


from equations (3.16) through (3.18) one obtains the 
following differential equations for f^(n) and 
H ± (n) (i = 1,2,3): 


1 i 1) ( ^i } + P \ = 5'i 1) ( r >) 
l( 2 ^ (5.) + L^^f.) = Pppn) 


(3.22a) 

(3.22b) 


The linear operators Ip ^ (i,j=l,2,3) are given by 

,3 


~ 

1 dn 




• n + b.f" 

o dn 10 


l(2) = JSL- + Prf „ - a^Prf ’ 


d7 


h 3 )= 

1 


'o dn 

i 2 

? 

_ A 

dn 


-2(1-Pr)(f --w + 2f g- + f^ *rr) + b^PrH f 

dn 


'*» d_ 

'o dn' 


(3.23) 


(3.24) 


(3.25) 
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w iinii f 3*^ — "b ^ 5 — kg — 3$ kii 0 fu.ncij'XOiis 

P^ 2 ^ are given by 


F'i 1} (n) = 4iP(H 0 -tf) + f”f 0 

for i=l,2 

(3.26) 

B'i 1) (”) = < 2 l- f oX f o- 2 V - % 

Hf'-(i + p )i[2 


+ | P(H 1 -2f;! , i) - Zp(H 0 -f; 2 ) 

for i=3 

(3.27) 

^• 2 ^(n) = PrH* f for i=l,2 

i v oc 7 


(3.28) 


lj 2) (i) = Pr^-U^X^-^) + PrSp?^) 


i*2 — -»*« 

+ 2(1-Pr) C ( f i ) + f i f i H for 1=3 (3.29) 


with d^ = 1.5 and dg = 2 . 


The boundary conditions can he written as 


f ( 0 ) = - 2 a 

0 

(3 ,30a) 

f‘( 0 ) = 0 
o' ' 

(3.30b) 

H^(0) = 0 (for insulated surface) 

(3.30c) 

or, H (0) = (for constant temperature surface) 

(3.30d) 

f’(n ) = 1 

o'' e 

(3.31a) 

H (n ) = 1 

o e 

(3.31b) 
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and for ? , EL (i = 1,2, 

f>) = 0 

3) 

(3.32a) 

f‘(0) = 0 


(3.32b) 

h’(0) = 0 

(for insulated surface) 

(3.32c ) 

or, 1.(0) = 0 

(for constant temperature 

surface ) 

(3.32d) 

ic e )= 0 


(3 .33a) 

V"e>= 0 


(3.33b) 


In order to determine the unknown constants 
Pi ,< "’i we insert the expressions (3.8), (3.10) and 

(3.11) in (2.50) to obtain the following expressions for 6: 


kill 

x 


(Po^ 



CVPiUrV 


3T72 

X 


+ { Pa( I 2- I o )+ &3 (I 2- I o )+ Pl 2 ( I 3- I l +I o ) Il k b 2 }/ x3 (3.34) 


where I = / 

0 

o 


CO 



< H o- f o ^ 

(3.35a) 

(H l- 8f o J l )dn 

(3.35b) 


(5 g -2f g )dn 


(3 .35c ) 
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for p . 


(3.15a) 


00 _ , „i 2 


On comparing (3.34) with (3.9) we get 


6 = Ini i 

0 ft>0 0 


6 1 = Pi ( :r - 1 ) 

0 


^2 ~ ^ 2^1 

o 


“ Pg(“ - 1 ) + P-j [I 


Os - 1 !) 


+ i3 


0 


(3.35d) 

( 3 ,36a) 

(3.36b) 
(3.36c) 
( 3.36d) 


Using (3.36) and (3.14) the following expressions 
are obtained: 


I (y -d C m^ii 3 1/2 i, 

8 P, 


4 v '~^' i- 2 

1/2 


P9 = 


(Y-1)(11I 0 -8I 1 ) 

8(3y+1) 

(Y 2 -l)[2Y(Y+l)3 i / 2 (13I o -10I ?; ,) 


p c (401 -401 ^+431 ) 
(52 I q -40I 2 ) 


(3.37a) 

(3.37b) 

(3.37c) 

(3.37d) 


The solution of ordinary differential equations 
(3.15b), (3.22a) and (3.22b) is obtained by the 



4:9 


shooting method of chapter IV. Except for the finite-difference 
approximations involved in the replacement of x-derivatives 
appearing in equations (3.1) and (3.2) whereby these partial 
differential equations are reduced to ordinary differential 
equations, the remaining procedure is similar to the one 
described there. 

3.4 Solution of the Steady Problem 

The series expansion solution discussed earlier has 
-3/2 

an error of order ^ \ and hence this solution is suitable 

only for high values of x. Since x can be taken as small 
as unity, the solutions obtained from the series expansion 
method could be in error for such values of x. In the 
present section a numerical method employing the complete 
tangent-wedge approximation of eq. (2.51) is presented. Since 
the complete tangent-wedge approximation is used, the obtained 
results are equally accurate for large and small values of "x, 

Eor the method presented here it is required to specify 
the initial values of f(x 1 ,n), H(x^,n) and p(x^) on an 
initial line x = x^. The values given by the series expansion 
solution on a line x = x^ (with x^ = 1/1024 corresponding to 

p 

X v 0(10 )) near the leading edge are used as the initial 
values. These initial values will be highly accurate due to the 
large value of x near the leading edge • 
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The initial conditions are, therefore 


f ,„ \ f / X P l^l^ n ' )k h Pg f g( ri ) + C ) 3 P 2 ( r i- )+ Pi P 3 ( 

fU^nj - f 0 (n; + —JJZ + — —- 3 — 


*1 


Xl 


(3.38a) 


H(x 1 , n) 


H f n ) P l H l^ n ^ k b P2^2^ n ' l+ S ) 3 H 2^ T1 ' >+P 1^3^ 11 ^ k b 
d o (n ' + —JJ2 + T — 


p( Xl ) = 1 + -pf + 


P l k b , V P 3 V 


(3.38b) 

(3.38c) 


X 1 Xl 


wnere 


x x 


The interaction between the boundary layer and the 
external layer is to be dealt with by an iterative procedure. 
Each iteration loop comprises of the two parts given below: 

(i) First the governing equations (3,1) and (3.2) are solved 
for a given p(x) distribution, in the region x 1 < x < 1, 

0 < n < n subject to the initial conditions (3.38a), (3.38b) 
and (3.38c) and the boundary conditions of section 3.2 and - 
(3.13b). 

(ii) Next a new p(x) distribution is obtained from the 
solutions f(x,n), H(x, n) of (i) making use of equations (2.50), 
(2.51), (3.5) and (3.38a), (3,38b) and (3.38c). 

The steps (i) and (ii) are repeated until p(x) 
converges within a specified tolerance. 



3.4a Solution of the Boundary layer Equations 


It can be observed from the series expansion solutions 
that the x-der ivatives of different flow variables, f ,H,p, etc. 
are large for small values of x and the x-dependence decreases 
in the down stream direction. Bor the same degree of accuracy 
we can use larger step sizes toward the trailing edge. Here 
in this work the following arrangement of x-stations with step 
sizes increasing in a geometric progression is employed : 

AX X = x 1 % 0 (10“ 3 ) (3.39a) 

AXg = 0 (10~ 2 J, Xg = + AXg (3.39b) 

Ax _j x j ”■ x ^ ^ ^ 3 3 ,4 , . . . , J (3.39c) 

Here k is a constant greater than 1 and J is an intege 
such that Xj = 1. 

In equations (3.1) and (3.2) the x-derivatives appear 

in the form 4x ~ . A finite difference formula for this term 

9x 

for the present arrangement of x-stations can be written as : 


4x 






(d D + V S 3-1 


+ E, £ o, j = 2,3,.. .,J (3.40) 

J J*~ <; 

where p. represents c(x.) and the coefficients D., E. are 
J J J J 

given by : 






j = 2 



52 


4(2k x +1) Xi 

nor — 


X 


3 

AX. 


3 - 3,4,,..,J 


(3.41a) 


Ej = 0,3 = 2 

4k 2 x j 
- x 3 

k +1 ax. * 

JL J 


3 — 3j4j»..jJ 


(3.41b) 


Expressions (3.41a) and (3.41b) correspond to two-point back- 
ward difference formula for j = 2 and three-point backward 
difference formula for j = 3,4,...,J. 

Substituting (3.40) in (3.1) and (3.2), we get 


p fV' + f 
* 3 3 3 


( 1+D . ) f . - (D .+E . )f . ,+E . f . „ 
3 3 3 3 3—1 3 3-2 


- f 3 [ ( VV f 3 ' ( W f 3-i + Vj-: 


and 


+ Q 3 H . = 0 


(3.42) 


P • 

*1 HV + H' 
P* 3 3 


(1+D .) f - - (D.+E . ) f . - 
v .V .1 v 3 3 n-1 


+ E 3 f 3-2. 


f'. 

3 


Dj H r ( Dj+Bj )H._ 1+E .E._ 2 


3 3- 


- (Pr-1 ) 

+ 2 Pj Tr“ 


f”) 2 +f f 
^ J J <3 


= 0 (3.43) 


In equations (3.42) and (3.43) the prime represents 
differentiation with respect to n and the function Q(x) is 
defined by 

2x d Pl 


Q(x) = § £l 


dx 


(3.44) 



are all known 


At any station x = x., f . H. 1 , f. 9 and H. 9 

functions of ^ and p., Q., D., E. are known constants. Thus 

u D <3 D 

equations (3.42) and (3.43) form a set of coupled nonlinear 

ordinary differential equations for f . and H. of the form : 

3 3 


i3 1 f m + p 2 ff" + ? 1 (n) f" + p 3 f 


• 2 


+ ^( r i)f* +p^_H = 


0 


(3.45) 


P 5 H" + p 6 fH' + 5 3 (n) H'+^f'H + ^(n) f' 


+ (3 


8 


(f ") 2 + ff m 1 =0 (3.46) 


where j3 1 s are known constants and uhe ?*s are known functions 
of n . The boundary conditions are : 


f k = 


4A. 

(f+4x fi } b = ~ 2a “ (j^n 


fp 0 M to p 


H, 


*1 ~ x 


Vx 


] 


H; = 0 or 


constant 


f'Oi e ) = H(n e } = 1 


(3.47a) 

(3.47b) 

(3.47c) 

(3.47d) 


By solving tne ordinary differential equation problem 
(3.45), (3.46), (3,47) at x = Xg, Xg,..,,Xj we can generate 
a numerical solution of the partial differential equations 
(3.1) and (3.2) in the region x, < x < x T , 0 < n < n . The method 
outlined here was first suggested by Hartree and womorsley [J59[[] 
and is known as the difference-differential or Clutter-Smith 
technique. in the literature |j5°H •- l^is has b .on us.jd extensively 
by Smith and Co-workers [50,61,6S] for the solution of 
boundary layer problems. Because 
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of the use of backward difference formulas for the x-derivatives, 
this implicit method is inherently stable jjpOJ . 

Equations (3.45), (3,46) and (3.47) can be considered 
as a fifth order nonlinear boundary value problem with three 
conditions specified at n = 0 and the other two specified 
at n = n g . Such a problem can only be solved iteratively. 

The most common method employed in the solution of such 
problems is the shooting method where the solution to tile 
boundary value problem is obtained by solving a sequence of 
initial value problems. The detailed procedure for doing this 
is described in Chapter IV. It may be mentioned here that 
the shooting method is an integral part of the Clutter-Smith 
technique jjSQ. 

3.5 Solution of the Interaction Equation 

After solving the boundary layer-equations the pressure 
distribution can be obtained using the equations (2.50) and 
(2.51) with 6 = 6(x) and p = p(x). First we introduce the 
normalized quantities 6(x) and I(x) : 

- - 

6(x) = M — — — T-jp- ■ x 4 (3.48a) 

* 0 xii 

00 r, 

I (x ) = / (H-f ' * ) dn/I (3.48b) 

o 

wnere 6 Q and I are defined m section 3.3. And from equation 
(2,50) p(x) can be taken as 
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p(x) = 


Prom (2.51) 


I(x) 

6 (x) 


1 r * ».T r r Y+ 1 


p( x ) = — i- r 

P 0 x i- 

which upon inversion gives 
K(x) = 


1 + YK { iji K + 


1 

(l+i K) 2 +l1 e 


(p P 0 x - 1) 


r 2 , - - „ N Y (Y + l) “V 2 

L Y + (p 0 p x " 1} 2™“J 

% 1/6 - X 1 / 2 

Jppg- 


or 


S(x) = 


where 


i/5 *-r-x 


1/4 


and 


S(x) = x 7 K(x) 

5 3 = p o * L * * 
How, from (3.5) 

K(x) = k, + M 


d6 

cTx 


and using relations (3,48a) and (3.52) we get 
S(x) = k b x 1 / 4 + \ (f 6 + x j|) 

- V 2 

% = 6 o X L * 


(3.49) 


\ 


(3.50) 


(3.51) 


(3.52 ) 


(3.53) 


where 
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or 


d6 = s ( x ) 
dx x % 


36 

4x 


k 


'b 

374 


x 


(3.54) 


when S(x) is substituted from (3,51), equation (3-,54) becomes 
a non-linear first-order ordinary differential equation in 6. 
This can be integrated as an initial value problem using fourth- 
order Runge-Kutta method. The initial value of 6 at x = 
is obtained from : 


OO W. 

f H(x 1 , n) - f ,<s (x 1 , n) 
6 = 6(x 1 ) = 


dn 


I o P^ x i^ 


(3.55) 


where f(x^, ), H(x^, ) and p(x^) are those given by series 
expansion solution at the initial line x = x^. Then having 
obtained S(x), relation (3.49) provides p(x) by using I(x) 
evaluated from the solutions of the governing equations (3.1) 
and (3.2). 

3.6 Results of the Steady Plow with Air Injection 

The results of the series expansion solution are 
presented in Tables 3,1 and 3.2 for Pr = 0.72 and y = 1.4, 

In the first table the various results are tabulated for 
several values of the injection/suction parameter, a, and a 
fixed value of the wall enthalpy, = 0.5. The second table 
contains results for a fixed value of a(= 0.3) and the wall 
enthalpy is given different constant values. Por comparison 
purposes the zero-order problem was solved separately for 



TABLE 3.1 


RESULTS PROM THE SERIES SOLUTION WITH INJECTIOE/SUCTIOI 



Pr = 

.72, y = 1.4 

, = .5 




a = .2 

a - .4 

a = . 5 

P 

II 

f 

* 

O 

a = -2. 0 

f 0 "(°) 

.4060 

.25188 

-.19849 

2.1587 

3.9971 

fl’(O) 

— 093 7 5 

-.07226 

-0. 16357 

-1.9843 

-3.9502 


-.17773 

-.11328 

-.09375 

-1.9423 

-3.898 

f 3 ( 0 ) 

0.11914 

0.03417 

.32226 

1.9328 

3.9296 

Hq ( 0 ) 

0.06927 

0.016911 

.002894 

.7 564 

1.4208 

H’(0) 

-.016301 

0.04026 

.03480 

-.6708 

-1.2811 

Hg(o; 

0.00835 

0.03764 

.036189 

-.65002 

-1.2635 

H£( 0 ) 

-.013256 

-.009733 

.007369 

.4 587 

.65495 

*0 

1.1754 

1.3252 

1.40498 

. 58447 

.38539 

I 1 

.36908 

0.3501 

.08751 

.47182 

.33776 


.41351 

0.6034 

-.2023 

.3882 

.34521 

J 3 

-.12489 

-.16632 

-.0782 

-.11594 

-.07040 

4 o 

.6954 

.73844 

.76034 

.4904 

.39822 

6 1 

-.92968 

-.8970 

-.93949 

-.6922 

-.6224 

& Z 

-.97217 

-.8133 

-.94259 

-1.510 

-1.1186 

6 3 

.2427 

,25029 

. 1243 03 

-.5452 

-1.6030 

Po 

.45706 

.5153 

.54632 

.2272 

.14985 

Pi 

1.3552 

1.2191 

1.001897 

3.5916 

5.0365 

Pg 

1.4997 

1.4934 

.82395 

4.4985 

10.729 

P3 

1.2682 

1.2058 

.6653 

1.4075 

.99818 
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TABLE 3.2 



RESULTS FROM 

TilE SERIES 

SOLUTION ’ 

tflTH INJECTION 



Pr = 

= .72, Y = 

1.4, oc = 0 

.3 



H b = 1.0 

H b = .75 

H b = .5 

H b = .25 

H b = 0.1 

f*(o) 

*4426 

.38023 

.31614 

,.25025 

.2096 

fj( 0) 

-.18359 

-.12817 

-.068359 

-.0062866 

.07031 

fl(0) 

-.2739 

-..20898 

-.13281 

-.022949 

-.000488 

^’(0) 

0.16601 

.-10937 

.05224 

.007049 

-.019287 

H'(0) 

-.06113 

-.007439 

.03953 

*078938 

.09842 

h;(oj 

.02173 

.019153 

.028287 

.032797 

.06308 

Hg(°) 

.02751 

.022769 

.031043 

.03359 

.03934 

HgCO) 

-.017475 

-.016695 

-.018896 

-;02361 

-.039399 

!o 

1.9224 

1.62917 

1.29879 

.92756 

.6783 

*1 

.6750 

.5474 

.45865 

.25123 

.07 672 

*2 

.7224 

.5925 

.31923 

.30157 

.08096 

I 3 

-.22859 

-.18767 

-.168995 

-.13361 

-. 18966 

6 o 

.8894 

.81875 

.73104 

.61779 

.5283 

6 1 

-.7125 

-.7805 

-.86591 

-1.0684 

-1.3303 

6 2 

-.5872 

-.69729 

-.92079 

-1.2475 

-1.838 

6 3 

.14636 

.17255 

« 16736 

.27521 

.22132 

p o 

.7475 

.6335 

.505029 

.36067 

.26377 

p i 

1.0981 

1.175535 

1.3386 

1.4672 

1.4999 

p 2 

.9408 

1.09587 

1.22087 

1.8486 

2.0872 

p 3 

.7893 

.92065 

1.00586 

1.4556 

1.8002 
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TABLE 3.3 

COMPARISON <71 TH THE SIMILAR SOLUTIONS OP COHEN AND 

RESHOIKO r63l 


p — • 5 f — * 6 y 

Pr = 1.0, a = 0. 

Ref . [63] j 

Present work 

! sJ = 0.2090 \ 

i b ! 

= 0.20774 

i f£ = 0.7946 | 

f£ = 0.7914 


(3 = 0.5, = 0, 6, Pr = 1.0, a = 0,0 and the results are 

tabulated in Table 3.3 against the similer solutions of Cohen 
and Reshotko Ml- There is quite good agreement between the 
results computed here and those of [63] • Table 3;4 shows the 
non- injection results by different methods used here against 
the values given by ref. [Ml. The variation between the 
results obtained in ref. [134] and those obtained here by 
Glut ler-Smith technique is due to the starting procedure 
adopted. In the work of Chattopadhy ay and Rodkiewicz Q>4] 
the solutions were started by fourth-order Runge-Kutta ’method 
and then switched over to the predict or- corrector technique 
for further integration across the boundary layer. Tne same 
step-size was used for starting as well as for further 
integration. Hence the accuracy of the results near the 
surface is questionable. In the present computations employing 
the shooting method a more accurate starting procedure 
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TABLE 3.4 

COMPARISON WITH THE NO-INJECTION RESULTS OE REE. 13d 


Pr = .72, y = 1.4, H b = .5 



Results 

Present Results 

Obtained by 



Prom Ref. 

m 

Clutter -Smith 
Technique 3d 

* Method of 

Einite Differences 

i?g 

fj(0) 


.6185 

.655514 


fj(0) 


-.3164 

-.33606 


f£(o) 


-.3554 

-.3597 


f"(°) 


.25585 

.26555 


H o<°) 


.1505 

.14985 


H|(0) 


-.06267 

-.05278 


Hg(0) 


-.04572 

-.04249 


fl'(o) 


.01615 

0.019339 


X o 

•ft -if 

.86725 

1.0314 

.86417 



.39369 

.50681 

.39215 


*2 

.37513 

.4618 

.37391 


X 3 

-.11225 

-.15649 

-.1117 


6 o 

0.59719 

.65146 

.59631 


6i 

-.99449 

-.8835 

-.99432 


6 2 

- 1.2628 

-1.0501 

-1.2652 


6 3 

0.30868 

.23246 

.31536 


Po 

0.33702 

.4010 

.33603 


Pi 

1.8134 

1.737 

1.8204 


P2 

2.2196 

1.9016 

2.2302 


P3 

1.8510 

1.5292 

1.87968 



*V al ae s obtained by Method ds s c r ibed in c Hap ter V 
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as suggested by Clutter and Smith [61] incorporating the 
Taylor's series expansion near the wall has been adopted. Table 
3.4 also contains the results obtained by a finite-difference 
method jl23] . Even though these values compare favourably with 
those of ref. jj/Tj "the present results obtained by the shooting 
method appear to be more accurate. The finite-difference 
method, however, appears better suited [23] in comparison to 
the Clutter-Smith technique d 61,65]] for problems involving 
surface mass transfer where skin- friction and heat transfer 
become small. 

In Pig. 3.1 and 3.2 the starting velocity and enthalpy 
profiles, f’(x^,n) and H(x^,n), obtained from the series 
expansion solution for x # = 8, = 20, H h = 0.5, y= 1.4, 

Pr = 0.72, 0^ = 2° and for different values of a are shown. 
These profiles are calculated at a distance of x^ = 1/1024 
from the leading edge which corresponds to the value of 
equal to 256. v/ith this value of X^, the profiles obtained 
from equations (3.38a) and (3.38b) would be highly accurate. 
Prom the starting station to the next one the step-size was 
taken to be approximately equal to 1/100 and the value of 
k = 1.4. This yields 14 stations from x = 0.0 to x = 1.0. 

The calculations in the n-direction were started with a step- 
size of 0.1/16 and this step-size was slowly increased to 0.1 
for use with the predictor-corrector formula! 

In Pig. 3.3 a comparison of the variations of p and 6 
along the plate obtained from the series expansion solution 
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and shooting method are given by Ivi = 20, x * = 8, H, = 0.5, 

Pr = 0,72, Y = 1.4, 9^ = 2° and for no-injection case. The 
agreement between the series solution and the shooting method 
of solution decreases along the length of the plate with the 
decreasing value of X because of the larger error terms in 
the series solution. 

Figures 3.4 through 3,7 show the variations of p,6,f-^' 
and H-^ along the plate for a fixed value of the injection 
parameter (a = 0,3) with different values of (=0.1, 0.25,0.5, 

0.75,1.0), M m = 20, x * = 8, Y= 1.4, Pr = 0.72, 0. = 2° and 

—6 

= 10“ . With the increase in the value of the wall enthalpy, 
the induced pressures are larger due to the greater displacement 
of the outer flow caused by the hotter fluid in the boundary 
la^er as shown in Figures 3.4 and 3,5, Fig. 3,6 shows the 
distribution of the wall shear function for different values 
of and the trend is similar to the one noted in Figures 3.4 
and 3,5. the wall enthalpy-gradient, however, decreases with 
increase in as shown in Fig. 3.7. This is due to the 
reason that the difference in enthalpy between the external 
layer and the plate surface exerts substantial influence on 
the wall enthalpy-gradient which is a measure of the heat 
conducted. 

Figures 3.8 through 3.11 present the effect of varying 
the injection/suction parameter, a, on the distributions of 
p, 6, f£ and along the plate for M m = 20, = 8, = 0.75, 
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Pr = 1,0, Y = 1.4, 0^ = 2° and. = 10**^. With injection (i.e. 
for positive values of a) p and 6 increase because of the 
presence of larger amount of fluid in the boundary layer. 

Because of the surface injection the velocity and enthalpy 
profiles are less fuller as compared to the no-injection case 
and accordingly their wall gradients will be smaller. When 
the fluid from the boundary layer is sucked out (i.e. for 
negative values of a) the boundary layer thins down giving 
smaller induced pressures due to the lesser displacement of the 
external flow as indicated by the lower values of 6. Because 
of the thinner boundary layer with suction, the velocity and 
enthalpy profiles are more fuller as compared to the no- 
injection case, resulting in very high values for their wall 
gradients. Hence whereas injection reduces the shear stress 
and heat transfer to the plate, these are increased by suction. 
However, the converse is true about their effect on the induced 
pressure. 

Figures 3.12 through 3.15 show the variations of p, 6, 
f£ and with 0^ along the plate of = 20, =8, Y = 1.4, 

Pr = 0.72, K 1 = 10~ 6 , = 0,5 and a fixed value of the injection 

parameter a(= 0.2). In Pig. 3.12 the increase in pressure with 
the increasing values of 0^ is due to greater inviscid compression 
with the larger wedge angles. 



Fig. 3.1 Starting velocity profiles* for different d. with %,= 256. 
M =20.0, %,* = 8.0, H h = .5,r=1.4,P r =.72,e h = 2 o .0 
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: '(x,, n.) 


Fig, 3*2 Starting Enthalpy profiles for different d with %j= 256.0 
Mqo =20, X^ = 8-,H b =.5,f^=.72,r=1.4 9 0 b =2.CP 
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Fig, 3 .6 Variation oi f b along the plate tor ditterent 
M oo = 2O,X L *=8^=.72,r' = 1.4,ot=.3,K ) = 1O»0 
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Fig . 3 .14 Variation of f b along the plate for different 0 b 
M oo = 20,>'X 1 *= 8, H b =.5 j oC= .2 ,T"= 1.4, FJ =.72,K 1 =10 







CHAPTER IV 


DETAILS OP THE CLUTTER-SMITH METHOD OE NUMERICAL SOLUTION 

4.1 Method of Solution 

The method of solution is similar to the one described 
in ref. QjfJ . The x-derivatives in the momentum and energy 
equations are replaced by backward finite differences, so 
that the partial differential equations are approximated by 
ordinary differential equations. Then the problem of solution 
is essentially to find the unknown boundary conditions at the 
wall that satisfy the known boundary conditions . This is 
done with the help of shooting method with a cut-and-try 
approach. This procedure is described in the subsequent 
sections. The momentum equation (3.1) and the energy equation 
(3.2) are interdependent (or coupled) and must be solved 
simultaneously . 

Of the several procedures possible for solving these 
equations simultaneously, the one used here is described below 
for a specified x-station. 

Initially an enthalpy distribution is assumed at the 
given x-station and then a solution to the momentum equation 
is obtained using the cut-and-try procedure. Values of the 
stream function f and the velocity function f ! from the 
first solution of the momentum equation are now used to obtain 
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a solution of the energy equation. This new enthalpy 
distribution is used to solve the momentum equation again. 

This iteration is continued until convergence of the solutions 
to the momentum equation is obtained to a specified accuracy. 
This procedure is depicted schematically in Pig. 4.1. 

The details of the method of solution of the momentum 
and energy equations at a particular x-station are given in 
sections 4.1 and’ 4.2 and the formulas required for performing 
the integrations are given in section 4.4. The method of 
integration is a predictor-corrector multi-step method that 
uses the Palkner multiple-integration extrapolation formulas 
and the Adams-type multiple-integration interpolation formulas. 

4.2 Solution of the Momentum Equation 

At any x-station on the plate the momentum equation (3.1) 
is first rewritten with the following substitutions: 

f = $ + ri 

f‘ = #* + 1 

(4.1) 

£ t 1 = $1 I 

f I I ' - $' ' ' 

Then the boundary conditions become 
at n = 0 

$ b = f b (4.2a) 


-1 


(4.2b) 
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| Assume f”(x,0) and a linear t 

' enthalpy distribution initially j 

T 

r - 1 

t 

Assume H'(x,0) or 

H ( x , 0 ) and solve 

Energy Equation 
■ with available solution 
of Momentum Equation 

j Solve Momentum Equation ^ 

Iteration ho. = 1 




Solve Momentum Equation j 


Solve Energy Equation 

I with the H(x, n)distr ib- i 


with the latest 

| ution obtained from 

1 J* H 

! 

solution of 

Momentum 

1 iteration 1 i 

: j 

j 

: i 

Equation 


Iteration Bo. = 2 

r — — « — 




Repeat iterations until 




< e 


where e is an accuracy input of 0(1CT 5 ) and v stands for the 
iteration number. 


fig. 4.1 flow diagram for solving the 
Boundary layer Equations at 
an x-station. 
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at n -* n 

e 

= 0 (4.3 

The substitution given by eq. (4.1) is made to reduce 
the round-off errors in the computer program. With the 
replacement of the x-derivatives by finite-difference schemes 
as explained in section 3.4a of chapter III, the momentum 
eq. (3.1) becomes a third-order non-linear ordinary differential 
equation. The boundary conditions for the equation at the 
surface and boundary layer edge are given by (4.2) and (4.3) 
respectively . 

Then the procedure is to solve (3.1) as an initial 
value problem with trial values of as a third boundary 

condition. Now the search for the correct value of 
is made by trying several values of until the solutions 

for are bounded in a specified region. Specifically, 

three solutions of the momentum equation are sought with 
trial values of Sq ' such that at n _ = ti is 

between the bounds of -Kg < ( n e ) < Kg. Both n and 

Kg are inputs to the computer program. Atleast one of the 
three solutions must be high ( $, ( r > ) > 0) and one low 
(4 ’ (n ) < 0). The Lagrangian three point interpolation 
utilizing these three solutions is then used to determine the 
correct solution that satisfies the outer boundary condition 
($'(n ) — 0 ) . The solution can be made as accurate as 
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desired by restricting the values of the bounds Kg. Typically 

K would have a value 1 for five place accuracy. 

& 

In considering the solution of (3.1) as an initial value 
problem, 1 is first determined from: 

= f n 4 ,,f an + (4.4) 

o 

Here the value of ’ is given by (3.1) and is the 

trial value. The integration in (4.4) is carried out by 
using the predictor-corrector technique at locations away 
from the wall. The Taylor’s series expansion with a shorter 
step-size is used to start the solution near the wall. The 
details of the integration procedure are given in ref. |_6l[ . 
Other quantities needed in the solution of (3.1) are obtained 
from; 



.1 , t ! 

* = J * 

0 

dn + 

(4.5) 


$ = / n «’ 

0 

dn + ^ 

(4.6) 

4.3 

Solution of the Energy Equation 



Similar to the case 

of momentum equation. 

here also we 


make the following substitution in the energy equation (3.2) 


H = G- + 1 
H’ = G-' 

H" - G" 


(4.7) 
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with this substitution the boundary conditions become: 
at n = 0 

= H b - 1 (4.8a) 

at n - 

G =0 (4.8b) 

With the values of the stream function and its derivatives 
known, the energy equation (3.2) is a linear second-order 
equation. The associated boundary conditions are given by 
(4.8a) and (4.8b). After replacing the x-derivatives by the 
finite-difference approximation, eq, (3.2) is converted into 
an initial value problem by the specification of a trial value 
of G b , Two trial values are used and the resulting solutions 
are linearly combined to give the particular solution satisfying 
the outer boundary condition. 

The integration formulas similar to those used for the 

solution of the momentum equation are employed to perform the 

necessary integrations. These formulas are used to obtain 

two solutions G^ and Gg of the equation (3.2). To avoid 

errors, one high and one low solutions are obtained such that 

G„(n ) > 0 and G 0 (n ) < 0. For better accuracies these 
p e e' 

solutions may again be bounded by: 

0 < G x (n e ) < K 3 (4.9a) 

° > Gg (n e ) > -Kg (4.9b) 

where K~ is an input to the program. Its typical value 
may be of order 100. 
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The two solutions G^ and Gg are now combined go give 
the correct particular solution G: 


where 


AG 


! + (1-A) Gg 


A = 


W 


W - 0 1 (n e ) 


(4.10) 


(4.11) 


4.4 Method of Integration 


This section is intended to give the details of the 
integration formulas employed in seeking the solutions to the 
governing equations. A special procedure involving the 
Taylor's series expansion is required to start the integration 
near the wall. The details of this are provided in Ref. [jol] . 
Here we consider the general situation (away from the wall) 
where the equations have been integrated up to r> n and we are 
interested to obtain values of $ and G and their derivatives 
at ( = n-n+ An). The extrapolation and interpolation 

formulas are used to approximate the integration indicated in 
(4.4), (4.5) and (4,6). The two-step procedure used in the 
integration is: 

(i) the extrapolation formulas are used first, with the 
known values of «’'* and ’ at the stations n, n-1, n-2 
and n-3, to obtain the values of $ , $ 1 , ' , and ' * at 

station n+1. The formulas employed are 
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^d b = 4’ + ■§ C55t"-59^:; + 37»-:' - 95 ;:- n £ 4 .i 2) 


4n+1) E = % + U 59^ + 37^1g - 9«;: 3 : (4.13) 


Wl) B = + «; + C3237; , -264^: i+ 159 7” 2 -38^: 3 3 


(4,14) 

where the subscript E denotes ’* extrapolated’ ’ . The 

5 

errors in these formulas are proportional to (An) or 
(An)^. Further details of these formulas are given in ref, (jplj , 
The value of b *’ 1 at n +1 can now be determined using the 
momentum equation and the extrapolated values of $* 1 , $ ! 

and $ . 


(ii) The interpolation formulas can now be used to 
determine more exact values of <5* * , , and $' 1 1 at the 

(n+l)th station. These formulas are 


t i 


n +1 


= ‘n 4 tKiiV + 


n+l'E 


19®i' '“5 


! » 1 

$ n~l + 


>’ * ' 
n -2 


(4.15) 


n +1 


n 


+ ill ra®'' + 19 4"-5®' * + e 1 ’ 1 

+ 24 Lr v n+1 + J n ° n-1 n-2 -1 


n +1 


(4.16) 


n +1 


^n + An 


+ 0 s84 


I I 


"n+ 1 + 17 1 0 n “ 3 6 °n~l + 7 $ n- 2 -Il 


(4.17) 

Finally the values of is obtained by again using the 

momentum equation and the interpolated values above. The 
errors in the interpolation formulas are much smaller in their 



87 


magnitude as compared to those in the extrapolation formulas. 
Also, tney are opposite in sign. 

The formulas required for the integration of the energy 
equation are similar to those given above. For details 
reference may be made to the work of Clutter and Smith 



CHAPTER V 


STEADY PLOW ON AN INCLINED PLAT PLATE '/TTH 
HETEROGENEOUS INJECTION 


5 .1 Introduction 

Equations (2.38) through (2.40) for the steady hypersonic 
flow on a flat plate inclined at a small angle 0^ to the 
free stream and with injection of a foreign gas may be written as 


5 n 


(x if| ) + f sf§ 

3n an 


- 4x 


af a 2 f _ af a 2 f 
3n 3x9n ” 5n 2 


(2 i || .1) 


H. 


■(— ) 2 
v 3n J 



l+c^( s-1 ) - 

i+c^(r^T7_ 


(5.1) 


P 


a_ 

an 


_x_ 

Pr 


|f +2(]?r-l) f£ 


2 

a f 

9 

an 


(Le-1 ) 






5C. 

an~ 


+ f 


8H 

8n 


-4x( 


a f aH aH af 
an dx ” an ax 


) = 0 (5.2) 


a *v dc . dc - pv 3c* 

£ 3n ^Pr Le + f an 4x %n 3x“ 


af ffi 

ax an 


) = 0 (5.3) 


Pur t her for a binary mixture (i=l,2): 


c 2 " 1-c l 


(5.4) 


Thus we need to retain the species concentration equation for 
one of the species only. 
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The equations (5.1) through (5.3) are solved by an 
extension of the Clutter-Smith technique [65] • The variations 
in the Clutter-Smith method for the case of heterogeneous 
injection are described in the next section. Other details 
are the same as those given in chapter IV. Before going to 
the method of solution we will consider some of the fluid 
properties required for the computations of the binary mixture 
of perfect gases in the boundary layer. The suffixes 1 and 2 
used for the fluid properties refer to the injectant (argonor 
helium) and the free steam fluid (air) respectively. 

The fluid properties required are c* , c* , c* } k*, 

^ P Pi Pg 

k*, k*, y*, y*, y* and D*g and some combinations of these 
that give the parameters Le and Pr . 


and c* of the mixture the following e xpressions 


are used; 


p* = 


K 


r M 1 M £ 


[ 2 > + M 2 


(5.5) 


where R* is the universal gas constant having a value of 

8.3196 x 10 3 — — 

sec^-K 



(5.7) 
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c* = 5.1988 xlO 3 — S- — (5.8) 

p He sec -K 


The variation of c* for air is small over the temperature 

ranae under consideration and constant mean value of c* ii 

^Air 

adopted to agree with the assumption of calorically perfect 
gas made earlier: 

2 

c* = 1.0383 xlO 3 — £L — (5.9) 

p Air sec -K 

The viscosity of the binary gas mixture is obtained from the 
following expression derived from Wilke's formula Q?6] : 


v> = 


-X- 

v l 


l+G 


x s 

12 T. 


+ 


11 * 

u 2 


l+G 


h 

21 T e 


where 


(1-c 2 )/M 1 

Cg/M2+( f _c 2 )/^i 


C 2/ M 2 

c 2 /M 2 +( 1 -c 2 )/ M i 


(5.10) 


(5.11a) 


(5.11b) 


G 


12' 



. 1/2 


M, 


v* 1/2 M 9 1/4 2 

{ 1+(~5f) (jp ) > 

V o 1 


(5.11c) 


Here u* 
spo cies . 


and ij* are the viscosities of the individual 
3 


M 2 


Air 




40 


29 
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The viscosities of the individual species are obtained 
from polynomials that fit values obtained from the kinetic 
theory. These polynomials jj>0] have the form 


yj = A+BT*+CT* 2 +DT* 3 +ET* 4 +BT* 5 (5.12) 

The coefficients appearing in expression (5.12) are shown in 
table 5.1. 


The thermal conductivities of the individual species 
are obtained from the Eucken relation : 


where 


c* . 
V3 


* 

k. 

3 


1 c* . 

i ( Q . „E.l _5 ) c * 

4 ^ o* 3 °vj 


/'“N 


is calcula ted by 


P 3 


= c* - 
V3 


R*/M, 


vz 


(5.13) 


(5. Id) 


The thermal conductivities of the mixture d-Stf-obtained 
by using the relation (5. Id) and replacing the individual 
species viscosities by conductivities. The binary diffusion 
coefficients D* g are obtained from polynomials that fit 
values of p*Dj 2 obtained from the kinetic theory. These 
polynomials j_20] have the form: 

p*D* = A+BT*+CT* 2 +DT* 3 +ET* 4 +ET* 5 (5.15) 

The various coefficients of the expression (5.-15) are 
provided in table 5.2, 
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Table 5.1 

Coefficients Appearing in the Polynomials for 
Viscosities of Individual Species, yV 

( 50K <_ T* < 4000 K) (Taken from Ref ) 


Coeffi- 

cients 



Gas 




He 


Ar 


Air 


A 

2.6346879 

10" 6 

-2.6511108 

10“ 7 

-9.25840 . 

10“ 8 

B 

6 .849132 

10-8 

9 .262746 

10~ 8 

7.698378 

10“ 8 

C 

-3.440510 

10“ 11 

-6.88502 

10“ 11 

-6 .078874 

IO - 11 

D 

1.129453 

10" 14 

3.898027 

10" 14 

3.274786 

10“ 14 

E 

-9. 0793378 

io- 19 

-1.038786 

10" 17 

—8 .290185 

10“ 18 

E 

-9.149432 

10"23 

1.036036 

10~ 21 

7.917685 

10”22 


Table 5 .2 

Coefficients Appearing in the Polynomials for Binary 
Diffusion Coefficient, in p*Dj 2 (50 K < T* < 4000 K) 

(Taken from ref. [68] ) 


Gas 


Coefficients 

He 


Ar 


A 

-4.341105 

10" 1 

-1 .445589 

10" 1 

B 

1.152878 

10“ 2 

3 .912013 

O 

l 

i C>1 

C 

5.751873 

10“ 5 

1.686750 

1 0 "5 

D 

-1.051026 

10~ 8 

-3 .399687 

10“ 19 

E 

2.1970975 

IO" 12 

6 .880437 

IO" 13 

E 

-1,983057 

IO" 16 

-5.965776 

10“ 17 
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5.2 Method of Solution 

At airy given x-station (denoted as m in the finite- 

difference scheme) along the plate, the fluid properties, c 1 

and G profiles computed from the previous (m-1) station are 

used to obtain a first solution of the momentum equation* 

The calculated stream function and its derivatives may be 

denoted by where the subscript NM represents 

the iteration number. Using the same fluid properties and 

the newly calculated functions the species and 

energy equations are solved. Since the fluid properties are 

known these equations are linear and their solution is 

comparatively simple. The computed profiles may be denoted 

by c- , cJ , G,, a _,, G T tqi?» where the subscript USE signifies 
TJSB J 'HSE iMbJi 

the iteration number. These solutions of the species and 
energy equations of USE = 1 are used to recalculate the fluid 
properties still using the stream functions and its derivatives 
calculated for NM = 1 earlier. The species and energy 
equations are solved again using these new fluid properties. 
These solutions, denoted by NSE = 2 are used to compute fluid 
properties once again. This process is continued until the 
solutions to the energy and species equations converge to the 
desired accuracy: 


where 


I g iJse “ 

6 p stxicl 


G NSE-ll < e 2 311(1 ^ C % SE ~ C 

are accuracy inputs, 
o 


| < e, (5.16) 
USE-1 ° 

Satisfaction of the 
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convergence condition (5.16) implies convergence of the fluid 
properties as well. The same stream function denoted by 
subscript NM =1 is used throughout this iteration loop* 

Upon satisfaction of the convergence condition (5.16), the 
momentum equation is solved a second time using the latest 
converged fluid properties and G and c 1 profiles. The 
computed stream function is denoted by NM = 2. This value 
of the stream function is now used again to solve the species 
and energy equations for obtaining the converged fluid properties 
as described above. This cycle of double iteration is continued 
until convergence of solutions of the momentum equation is 
obtained i.e., until 


m 


®NM-1 


< e. 


(5.17) 


To start the solution at the first station profiles of ^ 
and G are obtained from the series expansion solutions with 
air injection described in chapter III and the initial fluid 
properties for the iterative cycle are computed from these. 
During the course of computations of the heterogeneous 
injection case, it was found that four iterations of the species 
and energy equations satisfied conditions (5.16) if Gg and 
G 3 ranged from 10” 3 to lCT 4 . Consequently, this iteration 
loop was carried four cycles each time and the condition (5.16) 

was not built into the program. for the momentum equation 

-4 

the accuracy input was given a value 10 . 
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The flow diagram for solving the mutually coupled 
momentum, energy and the species concentration equation is 
given in Fig, 5,1, 

In general, convergence was substantially slower if the 
injected gas differed considerably (e.g, helium) in physical 
properties from the external fluid (air) as compared to the 
case when the injected gas did not differ much (e.g, argon) 
from the external fluid. Thus, whereas the Olutter-Smith 
method gave results with a fast rate of convergence for the 
case of argon injection, it failed to give the trial solutions 
within the bounds at the last station along the plate with 
helium injection. Hence an alternative numerical method was 
sought for in the case with helium injection which is described 
in section 5.5. Section 5.3 and 5.4 give outlines of the 
Clutter-Smith technique used for the case of argon injection. 

5.3 Solution of Momentum Equation 

v/ith the replacement of the x-derivatives by finite 
difference schemes as explained in section 3.4a of chapter III, 
the momentum equation (5.1) becomes a third order non-linear 
ordinary differential equation. Boundary conditions for the 
equation are specified at the surface and at the boundary 
layer edge by eqs. (2.56) and (2.58) respectively. The 
procedure to solve (5.1) as an initial value problem is 
the same as the one described in chapter IV. 
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In considering the solution of (5.1) as an initial value 
problem, ' is first determined from: 

M" = / (X0”)dri + X (5.18) 

Here X-^ is known and is the trial value. The value 

of ~~ (x$'' ) is given by (5.1). Similar to the procedure 
explained in chapter IV, the integration of eq, (5,18) is 
carried out using the predictor-corrector technique at 
locations away from the wall. The Taylor’s series expansion 
with shorter step-size is used to initiate the solution near 
the wall. Other quantities needed in the solution of (5.1) 
are given by 


4” 

_ Lhl'll 

X 

(5.19) 


J) 

= / f d n + 

(5.20) 


kJ 

0 


$ 

= J $ dn + 

t 

(5.21) 


o 


If may be noted that this particular formulation avoids 
the necessity of computing normal derivatives of X and thus 
eliminates a potential source of error. 

5.1 Solution of the Energy and Species Equations 

With the known values of stream function and fluid 
properties both energy and species equations (5.2) and (5.3) 
become second-order linear equation. Their solutions are 
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converted into initial value problems by the specification 
of (dG/on) b for the energy equation and c^ (and hence 
(dCp/Sn) b ) for the species equation, 'Two trial values are 
used in both cases and the resulting solutions are linearly 
combined to give the particular solution satisfying the outer 
boundary condition. 

For any x-station the equation (5.2) and (5.3) are 
rewritten to give: 


-($+n)G'+4x {($’ + 1) || - G* || } /p 


and 


3 n 


-(*+ n)cj+4x { (®'+l) 





( 5 . 22 ) 


(5.23) 


where 

n = n js'+2(Pr-l)(«’+l) «" + (Le-l) { j c 

p 3°i 

{(G+l) - («’ + !) > 

' l l ~ Pr Le an 


(5.24) 

(5.25) 


$ = f_n (5.26) 

G = H-l (5.27) 


With the substitutions given by equations (4.1) and (4.7) 
the boundary conditions (2.56) and (2,58) for the present case 


become 



98 


at n -* oo 

= 0 
G = 0 
= 0 

at n = 0 

«' = -1 


(5.28a) 
(5.28b) 
(5.28c ) 


( 5 . 29 a ) 


$ + 4x 


9J. 

3x 


b v b 


Hpe _ * ■ 

T V ( X ) J ^ 


Vo p 0 V x 


= (-^-)(l-e lh ) — — x (4 + 4x f|) 
Sn v u Le v lb' v dx' 


G. 


b 


= _°ib Cr - + x ] 


21 , 


(Y — l) 


£ 


(5.29b) 


(5.29c) 

(5.29d) 


In eq. (5.29b) the right hand side, in general, is a 
function of x. However, in order to have the results comparable 
with those obtained for the case of air (homogeneous) injection, 
el. (5.29b) is used in the form : 


$ + 4x ||r = -2a 
6x 


(5. 29e ) 


In the computer program, tile equations (5.22) and (5.23) are 
solved by programming a dummy equation of the form : 


• ($+ n) Q ! + 4x 


{(«'+!) |S 


Q 


i 

ax 


/p (5.30) 


where Z = Z(Q ! ,Q). 
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The associated boundary conditions are : 


at n = n 


at n = o 


Q = Q 


(5.31a) 


Q — Q 


'b 


(5.31b) 

Q' = (to be found) (5,31c) 

Solutions of the energy and species equations are then obtained 
by solving equation (5.30) with the conditions (5.31) where 
for the energy equation 


z = n (|| , G) 


’ Q = a 


Q, 


o 


Q b = G b 

and for the species equation 

- 

z = H (— ) 

Q = c. 


Q = 0 
^e 


C p, 


<$> b = (TT-eh^ 


x {$ + 4x f— ) > 

x b + v dx ; b 


(5.32) 

(5.33) 

(5.34) 

(5.35) 


(5.36) 

(5.37) 

(5.38) 

(5.39) 


It is to be noted that the specification of the wall 
concentration by relation (5.39) necessitates the recalculation 
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of aPber eacb solution ox the species equation. Hence G-^ 

varies with c* for a constant surface temperature. I*, c* , 
p b b Pb 

in turn, is dependent on tbe surface concentration, c^, 

obtained from the solution of the species equation. 

The method of solution of (5.39) is the similar to the 
one described in sec. 4.3 of Chapter IV. 

In the Clutter-Smith procedure employed hero and described 
in Chapters III and IV the x-derivatives are replaced by back- 
ward difference formulas and the partial differential equations 
tuns become converted into ordinary differential equations. 

Because of the use of backward difference formulas for the x- 
der ivati-ves, this implicit method is inherently stable. This 
method eventhough easy to apply and more efficient, has certain 
limitations especially, for the flow problems where skin-friction 
and heat, transfer become small, as in the case of boundary layer 
injection. Bor such c ases she derivative boundary conditions 

as n - °° may not be satisfied and accordingly the solution may 
not converge. With argon injection, however, the convergence 
was easily obtained. But when a very light gas like helium was 
injected, convergence was difficult even at the beginning stations 
along the plate and it did not converge at the last station. The 
trial solutions for the momentum equation could not be obtained 
within the bounds. The reason for this was observed to be due 
to the fact that when, at the last station, the value of ®"(0) 
was perturbed by an amount of 0(10 ) the value of $ at the 



At the start of a station, assume fluid properties from 
tne previous station. 
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Solve 

Species 

Equation 


USE = 1 


Repeat procedure in NM 
until NSE = NSI^ ax 


J? 


Solve Momentum Equation- 


NM 


Repeat NM-step until 


I - e 3 

Proceed to next station 


Eig. 5.1 Plow diagram for solving Boundary Layer 
Equations for a binary mixture at an 
x-station. 
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outer edge of the boundary layer was seen to change by an 

rz 

amount greater than 10 . Hence for the case of helium injection 
a finite-difference method based on Crank-B icolson implicit 
scneme was adopted, inis is described in tne next section. 

5.5 Finite-difference Method for Heterogeneous Injection 

To circumvent the difficulties which arose during the 
solution of the boundary-layer equation with helium injection 
a finite-differences method described in [25] is briefly presented 
here. In this method, the governing partial differential equations 
are transformed into new coordinates with finite ranges by means 
of a tr. msformation whicn maps an infinite interval into a 
finite one. 'Ihe resulting equations are solved by converting 
them into a matrix equation through the application of implicit 
finite-difference formulas. Ihe matrix equations have a 
tridiagonal character and can be solved on a computer by making 
use of a suitable algorithm [69,23,70]. Only those details 
which are pertinent to the present problem are given here. For 
further details a reference may be made to the original paper 
[23] . The distinct advantage of the present method over the 
Olutter-Smith technique is that values of the dependent variables, 
rather than their slopes, are employed to obtain solutions. The 
associated disadvantage with this method is that there might be 
slight oscillations present in the quantities like 9 f/9n and 
aH/dri etc, due to their computations from the finite-difference 
approximations. 
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First of all, the equations (5,1) to (5,3) are transformed 
from the (x- ^-coordinates system to the (x- s) system thereby 
replacing the infinite interval in the n-direction into a finite 
interval, The transformation employed is : 


S 



(5.40) 


For a system of finite-difference equations with a fixed number 
of S nodal points and a fixed interval AS, a ^ is used as a 
scaling factor to provide an optimum distribution of nodal 
points across the boundary layer. A certain computer experimen- 
tation is done to achieve tne optimum results. The following 
transformation equations are obtained from equation (5,40) 


_ 95 a 5 _ ( 1 r'l 

on as 6n a l^ 1 ~ v 



and 




a l 


Z 


a? 

as 


(5.41) 


(5.42) 


where Z = a^(l-5) 


(5.43) 


Using these transformation equations and using F = f ’ , equation 
(5.1) through (5.3) become 


p Z 2|||| + p XZ ( Z ||-a 


8F 

i as 


) = -f z 


3F 

as 


. af aF v 
4x < z ax W ~ F 


aP ) + p(% 


ax 


p 


dx 


o r- 1+ e -, ( s- 1 ) ■ 

DtH-sn — L —— 

*— 1+c ^(r— 1 ) 1 
(5.44) 
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5 z * [ft <j& 


P 3c 1 1 

(H-P ) 


sS. x i \t? 3-^ ' ■ 1 

05 


+ 2(Pr-l)P H + (Le-1 ) { 


+ P (p^) 


Z 2 LL _ a z M 

55 2 1 3C 


(r — 1 )c-j+l 
QH 


> 


+ f + «S(M-1> fZ 2 (||) £ 


+ 


^ 2 Q 2 ^ 3F \ x , ^2 9Le r-1 

-— - - a x Z F > + Z ~~~ { — “ — rr~r > 


35‘ 


35 c 1 (r-l)+l 


o 3 c -i 

(H-P 2 ) + (le- 1 ) { 


, 8 „ 

r-i p p o c. 

— — > (H-F 2 ) (Z~ — 7T 

c 1 (r-l)+l d 5^ 


5c 


5c. 


- a 


Z Zl, + (le-1) (— Z 2 (ff -8P m 


1 " 55 


c^r-D+l 


35 


/T , Wttto 2 w 2 r "(r-1) , / 8c l^ 2 1 

+ (Le-1)(H~P )Z { — ?} ( 77 -) 


( (r -1 )c 1 +l ) 2 35 


. f 7 aH, , 7 m af x 

= -f Z ^ + 4x (F " Z 35 


(5.45) 


5 5 v Pr ' 5 5 
Oc 


Pr 

5c. 


- -z f ^ + 4x (I ri _ 2 M £l) 


3 5‘ 

01 
3x 5 5 


- 7 2 a / XLe ^ Qc i , = xLe / 7 2 3 " C 1 „ z dC l ) 
P L TT 77~ + P “Pr~ ( Z “ “ a l Z 5 5 ' 


where 


f (x, 5 ) = / 7 h-5 + f (x, 0 ) 


(5.46) 


(5.47) 


l he boundary conditions become 


at 


0 
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F 

= 0 

H b< : 

z) 

■[ 

r c lb 

f + 

4x 

Of 

Ox 

= -2a 

0 c^ 


(1-c 

lb ^ 

” 


a 

1 


w] 


(y-i) K 


(5.48a) 

(5.48b) 

(5.48c) 


v yLe'b Pv, 


(f+4x (5.48d) 


and. at 5 = 1 


F = 1 
H = 1 


0^=0 


( 5 .49a ) 
(5.49b) 
(5.49c) 


5.6 Finite-difference Equations For the Hen-Similar Problem 

The mesh-point diagram for the Crank-Hicholson scheme 
is given in Fig. 5.2. tfitn the N nodal points in the 5 -dir ection, 
we have 
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(H - 1) Ag= 1 (5.50) 

She following finite-difference equations are used to 
reduce the non-similar boundary layer equations to a system 
of linear algebraic equations 


G = G 


G = 
x 


e 

s e - t 

AX 


G 5 " iT? ^ G e “ G a + G f ~ G d^ 


G 


55 


(v so . b ♦ Ga + Sf 


M - Mr 

G1'.I = Gr - 

x b ax 


(5.50a) 


(5.50b) 


(5.50c) 


2G e + G d^ ( 5 « 50d ) 


(5. 50e ) 


G 5 M 5 


8( A ?)' 


(G e -G a H M f- M d ) + ^V G dHM e -M a ) 

( 5 . 50f ) 


G. 


4 (a zy 


(®e - G a> ( G f " %) 


(5.50g) 


G ~ = G S xG e 


R« 


R 


R n - R 
c a 

2(A?)“ 


R c “ 2R b + R a 


55 


(a ? ) ( 


(5. 50h) 
(5.50c) 

( 5 . 503 ) 


Equations (5.50b - d) are used to approximate the unknowns 
and their derivatives* These are the Crank-13 ic ho Ison formulas. 
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The equations (5«50e - h) are used to make the unknowns to 
appear linearly. The central-difference equations (5.50i - j ) 
are for the derivatives of fluid properties. 

5.7 Method of Solution 

To obtain a numerical solution, equations (5.44) through 
(5.46) are approximated by tne finite-difference equations 
described in the previous section. The resulting equations are : 


1 


d 


_ 2 — 2 — 

pZ Xg pXZ pXa^Z zf xZ 

4 ( A 5) 2(AC) 2 4(a$) 4AC ( AC) 
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r p x z‘ 


(A5)‘ 


E b tpl f E 1 1+Cl - (r _i) S J 




t-p Z 2 X f p X Z 2 pXa-,Z zf x Z f -j 
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4 A 5 4 AC 
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Zf 
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x Z f 


4 A C 


a' 4 A? <*c " 


— £ (F -P ) 
A? c a 


— F- 

Ax b 


1 + c.r (s-1) 

F- + p(~^ ^ - 1) H r { — > (5.51) 


-2x 

dx 


1 + c ib 


1 - X Z 


H, [-p Z 2 (j^) 

p (X/pr ) Z^(Le-l) (r-1) 
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2 

Z^P 
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£LZlI 5 l 
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+ 2f + x2f x ~ 
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Ill 


In equations (5.51—53) following the first iteration 
quantities not subscripted with V are evaluated at 
(n , m + 1/2) by averaging quantities at (n,m) with those 
at (n,m+-l). 


how for the species equation both 0 .( 0 ) and 0^(0) 
are unknowns at the wall. A special formulation is required 
to express c^(0) at the wall. The procedure is to express 
c^(0) in terms of c^(A£ ) as outlined below; 


Prom the boundary condition equation (5.48b) we may write; 


110 ,( 0 ) 1 . 


~ [i-cho)]^) 

<*1 U 1 J Le b ( UP )^ 


(5.54) 


The species equation (5.46) specialised at the wall gives 


p a 
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~/XLe \ 
v Pr ) 


Co (0)3 + p(i^) 

b 1 s 


Pr 'b a l 


JCi(°)^ 5 ~{c 1 (0)}^ 


= -a x f b C ° 1 (0)l^-4xa 1 Z£(0)1 X [c^O)^ (5.55) 

How, from the lay lor* s series expansion we have 

2 

CjUe) = opo) + 4 5 Co 1 (o)3 e + Bi<o) 3 k+ ••• 

(5.56) 

Substituting for ^c^O)] from e< l* (5.54) and for £<^(0)]^ 
from eq. (5.55) in the above equation, we obtain 

c 1 (A?) A s +B s 

°i (0) - + 


(5.57) 
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vhere 
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( 5. 58a) 
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a 1 v Le ; b 


0 Pq x 
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r / XLe 
u Pr 



P(p r 'To 



/pCH 2 -)] 


a 


2 

1 


(5.58b) 


5,8 Starting Profiles 

Phe solution at the first station on the plate is started 
by using the series expansion solution developed for air injection 
in chapter III. To use these solutions, only the air infection 
is used at the first two stations in the x-direction i.e. upto 
x=0,02. Beyond this value, the foreign gas injection is 
employed. It may again be mentioned here that initiating the 
eolations at x=0 as done in ref. D=fj and other worhs is not 
correct, especially, for a high speed problem as pointed oat 
in section 1.2 of chapter I. An alternate method could have 
been to develop a series solution for the binary gas-mixture. 

With the variable fluid properties for this case, this is 
somewhat more complex task. 

In this section we recast the equations of the series 
expansion method for air injection in terms of the variables 
The resulting solution is then utilized for 


x and £ . 
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specifying the initial conditions in x for obtaining the 
:inite difference solution described in the previous section. 

Through the use of transformation equation (5.41) and 
[5.42) we get the following equations for (3.15) and (3.22). 


? 9P_ 9? o 
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In equations (5.59-62) tJae following substitutions are also 


made 
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, -« _ _» _ _f 

F o = f o» ^1 = f l» -^2 = f 2 ^3 = f 3 » where primes 

denote differentiation with respect to r\ . The equations 

(5.59—62) correspond to similar form and a special scheme for 
finite differencing is to be adopted as given below. 
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fig. 5.3. Mesh-point diagram for similar solutions 
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where superscript v is the iteration number. 

Substituting (5.63) into equations (5.59—62), we obtain 
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In equations (5.64) through (5.67) 'the quantities not 

T _ ¥ 

subscripted with b are taken from the previous iteration. 

And after each iteration the functions f Q and ft (i=l,2,3) 
are updated using integral relations similar to that of eq. (5.47). 


5.9 Similarity Form 

The similarity forms of the equations (5.38) through 
(5.40) are solved for comparison with the existing results [if] . 
These equations are obtained by taking the value of p as 
zero and p as unity: 
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Substituting from the finite difference equations 
(5.63) into the above equations we obtain: 
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For solving eq. (5.73), the expression for c^(O) 

(i.e* equation (5.57)) is to be modified by removing the term 
involving x-derivative. 


The similarity forms of the equations for homogeneous 
injection case are also solved for comparison purposes [ll] . 
Applying the transformations as before, the similarity forms of 
equations (3.1) and (3.2) with §=0 yield: 
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The f inite—dif f erence form of the above equations may be 
obtained by using the expressions (5.59); 
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5.10 Results of the Steady Heterogeneous Injection Problem 

Computations were carried out for a sharp wedge of semi- 
apex angle 0^ =2° at a free stream Mach number of 15. I he 
value of used was based on a constant surface temperature, 

j of 1948 .5K. Corresponding to an altitude of 30480m 
and standard atmospheric conditions, the following values of 
the various free stream variables have been used in the 
numerical computations: 


I ,c 

CO 

= 216.5 

K 

£ 

= 295.0 

m/s 


= 4440 

m/s 


In figure 5.4 the variations of p and 6 along the 
plate for a constant injection value of a=0,l for air and 
argon is given with 1^=15, = 8, y = 1.4, 6^ = 2°, and 

= 9 . Since air is lighter than argon, the induced pressure 
and displacement thickness for argon will be lower than those 
for air injection with the same value of the injection parameter 
a as shown in this figure. 

Figures 5.5 through 5.8. show the variations of p, 6, 
f” and for different values of a (for argon injection) 

against the corresponding no-injection curves for 15, 

X L * = 8, = 9 , Y = 1.4, and 9^ = 2. Similar to the case for 

air injection, the induced pressures and the displacement 
thickness are larger whereas the wall shear function and the 
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enthalpy gradient are smaller with injection. 

Fig. 5.9 shows the variation of wall enthalpy with 
argon injection against a non- in jection case. The wall 
enthalpy decreases as a increases because the value of 
c* for argon is less than that of air. Fig. 5,10 gives the 

Jr 

variation of argon concentration through the boundary layer 
at x = 0.5 for different values of a (=0.1, 0.15) with 11^=15, 

= 8 ’ T b = = 2°» Y = 1.4. Figure 5.11 gives the wall 

concentration of the injectant for different values of injection 
(a = 0.1 and 0.15) along the plate for = 15, = 8, 

= 9 , Y =1.4 and 6^ = 2°. In obtaining the results of 
Figures 5.4 through 5.11 the Clutter-Smith technique employing 
the shooting method was adopted. 

Figures 5.12 and 5.13 contain the distributions of p , 6, 
f^’, and along the plate for argon injection by two 

different methods: (a) Clutter-Smith (shooting) method (b) finite- 
differences method with = 15, x^* = 8 j ^ = 9 5 = 2°, 

Y = 1.4 and a = 0.1. As pointed out earlier, the Clutter-Smith 
technique did not work for the case involving the injection 
of a light gas like helium. The finite-difference method 
employed here to overcome this difficulty does give quite 
accurate results for both heavier and lighter than .air injectants. 
This is evident from the comparison of this method with the 
Clutter-Smith technique for the case of argon injection as 
contained in Figs. 5.12 and 5.13. The accuracy of the finite— 
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difference method for the cases with air and helium injections 

f 

is shown in figures 5.14 through 5,16, These figures contain 
comparisons with the existing results [11,15} . As can be 
easily noticed, there is very good agreement between the present 
results and those of refs, [if] and [ 16 ] in all these 
figures. The results with c arb on- di- oxide injection in fig. 5.12 
have boon included for comparison against the argon injection 
as both gases have comparable molecular weights. It may be 
noticed in this figure tuat due to the difference in variables 
employed here and in ref. [16}, proper scaling has been done 
for the purpose of oomparison. 

In figures 5.17 through 5.2 0 are given the distributions 

of p , 6 , f ^ ' and for different injectants by the finite- 

differences method with M m = 15, =8, 9^ - 2°,Y = 1.4, 

♦ 

= 9 and a = 0.15. The pressure and displacement thickness 
increase according to the decrease in molecular weight of the 
injectants as shown in figs. 5.17 and 5.18, This is in 
accordance with the observations of figure 5.4. The f^' 
curve for helium in figure 5.19 shows a different behaviour. 

For the value of the injection parameter, a, considered here, 
there is an increase in the shear function as compared to the 
air and argon injection cases. This behaviour, however, can 
be explained if a reference is made to fig. 5.15 described 
earlier. It may be noticed from this figure that at first there 
An an increase in the value of wall shear function with the 
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increase in infection parameter, a, up to a value of cx = 0.1625. 
Beyond, this tne shear function decreases and eventually its 
value becomes lower than the case witn air and argon injections 
for a > 0.2l. The results for the enthalpy gradient at the 
wall are contained in Fig. 5.20. Tne trend of the results is 
similar to the one shown in Fig. 5.19. 

Tne present analysis uses values of a in the range 
-2 < a < 0.5 with negative values of a indicating suction. 

This is approximately the range covered by Zien J~ 14= . However, 

the detailed results have been obtained for a value of a = 0.1 5. 
This corresponds to a value of f^ = -0.3 shown in Fig. 5.15 
wnere a comparison with the Baron's results Q 15^ nas been made- 
The present analysis loses its validity in the vicinity of 
separation point where f^' approaches zero. Figure 5.15 shows 
the range of a for different gases before the separation occurs. 
The present analysis was limited to the values of a well before 
the separation point is reached. 

Finally, Figs. 5.19 and 5.2 0 are not the correct indicators 
of wall shear and heat transfer as is the case in absence of 
surface injection. Defining the shin-friction coefficient as 

°£ - ^*/! ■>: < 5 - 76) 
one may write, in terms of the transformations of Chapter II 
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Or , 

Therefore, 


c f * x ~ 3 / 4 p(f M -k) 


c f “ c f 


<p £ i% 


- 1 


(5.78) 


(5.79) 


where the subscript denotes the values without infection. 
Similarly, by defining the dimensionless heat transfer rate as 


= - q w 
^-w 1 ,, *3 


(5.80) 
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we may write for the Stanton number as 
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(5.81) 


or, in terms of the transformed variables of Chapter II, 
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(5.82) 


Figures 5.21 through 5.24 give the percentage changes in 
induced pressure, skin-friction coefficient, Stanton number, and 
boundary layer thickness with injection of various gases. These 
figures suggest that the lighter gases are more efficient 
coolants. However, they induce larger pressures. 
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Fig. 5. 8 ARGON INJECTION-Shooting method (Clutter- Smith) 
Variation of along the plate for different d 
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Fig. 5.11 ARGON INJECTION -Shooting method (Clutter- Smith) 
Variation of C 1t? along the plate for different d 
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Fig. 5. 15 Variation of f b with f b for different injectants i 
solution with M oo =4.0,T*=1647°R , T b = 1167°R 

Te = 392°R,u (2 = 3880 / /sec 

O BARON [15] 

PRESENT WORK (FINITE- DIFFERENCE METHOD) 
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CHAPTER VI 


UTILIZATION OF RESULTS OBTAINED WITH SURFACE 
MASS TRANSFER 

6.1 Calculation of Important Parameters 

In Fig. 6.1 the various forces acting on a slender wedge- 
wing of semi-apex angle at an angle of attach 5 is 

shown. The subscripts T and B represent the top and bottor 
surfaces respectively* 



Fig. 6.1. Aerodynamic Forces Acting on a 'Wedge-Wing 

The total normal force permit width acting on the two 
sides are 
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J (pj^g) dx * 
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( 6 . 2 ) 
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The tangential force per unit width on the two sides 

are given by 


( f Pt,b 


o 


L* 

f (?*) T> b dx* 


(6.3) 


where t* is the shear stress on the surface given by 


x* = (u* -d „ 

( ay^V =yj 

Using the variables of chapter II, we may write: 
T * _ 1 *2 ”^ p o / r x3/2 - -3/4 .it i 

T - 2 P o° u - JJST < X L*) P x 7 1 Mn =0 


(6.4) 


(6.5) 


Hence , 


3/2 

,B _ V"P- 1 

| p* u* 2 

2 oo OO 
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Now, the lift and drag coefficients can be obtained from 
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where ©g - ©>, + a 
and e T = e b “ “ 

Using the expressions (6.2) and (6.6) in eqs. (6.7) and (6.8) 
we obtain: 
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where S-^ and Sg represent the integrals 
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The aerodynamic centre is located by 
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In the non-dimensional form 
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(6.15) 


The heat transfer rate to the wedge per unit area at 
any point of the surface is given by 


3 rp 3 C js | 

k* + (c*. - c* )T* p* B* ~ 

pi p2 12 3y*_j y=y£ 


(6.16) 
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Using the transformation of chapter II, we obtain: 
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If y* denote the total heat transfer rate per unit width 
to any one side of the wedge 
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In the case of air injection the integrals S g and 
S 4 become 

<v ai r = {Y*- 3/4 ? £ "3b to (6 - 20) 

(» 4 W ~h{ x ' 3/4 f } * dx (6-21) 

The integrals S 4 (i=l,2,3,4) appearing in equations 
(6.H), (6.12), (6.15) and (6.19) have been evaluated by the 
Simpson’s rule for unequal intervals as described in appendix B 
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I* or the evaluation of C^, etc, the flow governing 
equations arc solved twice for each value of the angle of 
attach: once for the top surface chid second time for the 
bottom surface of the wedge. Obtaining of these results from 
the numerical integration of the governing equations, therefore, 
requires substantial computer time. The series expansion 
solution cai also be used [?l] to obtain these results without 
no. ding much computational time. She utility of these results, 
however, is limited to large values of the interaction parameter 
X » Appendix 0 contains the expressions for various aerodynamic 
coei'f ioionts obtained by utilising the series expansion solutions 
for the case with and without air injection. 

6 • 12 Results of the Aerodynamic Characteristics with 

Surface Injection 

Pig. 6.2 shows the comparison of the values of °I> °J> 
and JLi/j) obtained by the Clutter-Smith method and by using 
the series expansion solutions for various values of the angle 
of attack for the case of air injection with = 20, 8» 

h = 0.5, i J r = .72, y = 1.4, 0 b = 2° and a = 0.5. She 
agreement between the results obtained from these two methods 
appears to be good at small values of the angle of attach, ai 

Table 6.1 gives the values of the integrals S^i^^^jd) 

and C, , C n , L/D, x, „ and Q (Q is the total heat transfer 
Jj 7 D 7 etc 

rate to the wedge por unit width obtained by adding the values 
for top and bottom surfaces of the wedge given by eq. (6.18)) 
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for nil ' ngle of attack of 1° for both homogeneous (air) and 
heterogeneous (helium and argon) injection cases with = 15, 
= = Y = 1.4, a = 0.15. The calculations of the 

values of C^ > 0^ and L/f for various values of the angle of 
attack was not performed due to the prohibitive computer time • 


fable 6 .1 

Values of the Parameters for the Calculation of tho 
Aerodynamic Characteristics of a wedge-wing with 
Surface Mass fransfer (11^=15, x L * = 8, f^ = 9,y=1.4, 

a=.l5, 9 b =3°, Q fj} = 2°,0 d = 4° and Pr = .72) 

Angle of attack, a=l° 


p ar aiue ter 


Injectant 


Ar 

He 

Air 

S lf 

3.18698 

3.67411 

3.23367 

S 21 

2.20448 

2.01544 

1 -.76914 

S 3f 

1.21723 

1.42376 

1.24557 

S 1B 

4.02433 

4.46287 

4.22231 

S 2B 

2 .30056 

2.11976 

2.00430 

S 3B 

1.61833 

1.80220 

1.72011 

^4B 

1.42956 

1 .20748 

0.97002 

°L 

0.01004875 

0.009469498 

O'. 01193239 

°JD 

L/D 

0.01971726 

0.01908277 

0.01749107 

0.5096424 

0.4962328 

0.6821985 

nC 

0.4921834 

0.4927840 

0.4903794 

Q 

0.00450726 

0.0020360 

0.004127421 
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CHAPI’ER VII 

» 

UlS'IEADI CASE OE HOMOC-E1EOUS BlJECflUN 

Cni;j chapter wo propose to consiier toe problem of 
'■ U u.t -1 ti tent of truj bean imy layer when the inclin- tion of 
the o J. J rj Lrueulsivoly changed from the value 6, . to b v 

Vhu [‘ioo'dLa.j followed is -inscribed h- ; re under. 

0 t Op j- * 

C'loul. to the bo un lory layer cn or acteristics — f(.x, ) and 
:i(x, ) an l the surface prussuie; distribution p(xj separately 
eorr espon im , to tn 0 incline cion u ;nd 6, let 

D j 1 U y 1 

tuju be : coign, tod as id, p i? f f , Up and . Also 

c 'deal, at „ too surface* pressure listr lbution with K of 
equ Lioii (b.bj c; Iculatc 1 with k, Q corresponding to u. ^ 
ani o oorr espon ling to 9^ Call it Pj_f* 

Et- . o :i - 

Ait: Cue unsto s / , -o verging oqu* tions ^n l use the 
ini fix J *;ou 1 it i.Ui.uj 

“ 1 i 

-^ x ’ r '^=c “ lf i 
P( x ^=e = p i-f 

Stop b . 

Solve tun - unsteady ,.overnin te equations 'aid calculate 
tno turuj it'puirol for f,ii and p to attain the values 

1 ,p y • 'Xi L P -| * * 



in otiitiL* worar.:, it is propose d, to calculate the time 
r<:iUi.iO'i ior tin- boun d," X'j is yt-o? to change from tne steady 
r.tate correspond in.'.; to plutw inclination 0.^ .. to that 

_ *- — “ — - D , 1 — — 

cou'i.i'vmi ;iir, to bin* plate inciinet ion + , if trie pressure 

’’" Tn " 1 " u J rT ' r im 11 '‘""“""TI" 1 " '** "**” """ 1 "***"'“*' ~ r nr - "" l ' - -mn-mtmmmmmm [} J i. wr. 1 . 111 . unr i rii-Tiirni'n-- ’ r, "’ 1 " r ” ' 

dip trim dm 0 1 - t °n tn , plate is impulsivel y c hanged from the 
one cor.. - -:n,eu. iin . to o,. , tc that coir esponning to 0^ ^ 

(with h corf o spend in, to h ■). 

\ fJ ? i' 

7.1 r'l.cvi io'/.,rnmd Equations and tnc Initial Conditions 

i *n>. ovornmg o-. nations for this ur.se may be obtained 


from .l:i j - nd ( 2 .3 -j ) of Un.-ipter rr by putting c^ ■=• 0 

in cur. so u-u tions : 


f 1 


;* U 


— *i.X 



r < « _ 



+ tj < 12 iE. 
P ^ 

, 4x afl 

+ 4x at 


-1; (E-f 1 ) 

(7.1) 


1 r 


i'» + i n 1 


_ l , djj 


* oi' 


g-; = II (W-J U " 2 + f'f") 
5il 


u :: 


+ 'lit 


et 


-- §1 (-«'! 


(7.2) 


•« rn.r p ( :: , z ) i j .,1 v on by expr eosion (2.51;. 

fno steady state, solution of o nap tor ril» corresponding 
to u = u , decor im the flow on the plate prevailing at 

j) b y 1 

t < u. it tin: inst uit of the impulsive change (t=0), the 
spatial diotr ibution of th- temper' tore and velocity field 
w j.tn ry; ,'i.ct to tno plate is supplied oy this steady state 


;ne ch;in e in the inclination of the body 


solution. hut 

ctmngosj tn- effective shape of the body which in turn changes 
tnu pressure distribution. 
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Tno initial conditions in x are specified at x = for 
^ 2 l Pn ^lio steady problem, the series expansion solution 

I or is assumed to give the solution at x = for 

i 2. ( ~ >m Itiis assumes that the final steady state at x = x^, is 
reached instantaneously* Because of the close proximity to 
tiio leading edge, the flow at x = x^ is extremely insensitive 
to the inclination and the initial and final steady state 
solutions differ by less than 1 percent, fhus the initial 
conditions in x are given by 


(i) li'or x = x 1 at t > 0 

4.\ _ * / \ . P 1 k b,f 

n? ^ ) — £q ( n ) + r n ' 1 ,r 

Xl 


+ 


■ p 2 h^-1 1 t p i 


(7.3) 


H 1 (n) p ± k, f 
ii(x 1 ,n,t) = II Q ( n) + 


p 2 Hg(n) + Cp 3 Hg(ri) + Hg(n) f 


, v - . p l k b 7 f P 2 + p 3 
p(x 1 ,t) = 1 + — 337 J- + — 


(7.4) 

(7.5) 


H 


where x*. = X * /]fx 1 and k-u £ = M m 0^ f 

X jj X 3 3 

(ii) For x > x ± at t = 0 

f (x, n , 0) = f ^(x, n) 


(7.6) 



J-CKt 


ii( x ;*i, u ) = i/ U ' ) (x, v u 


wtit-r*. 


(7.7) 

^ and ^ ^.ct: the steady state solution corres- 


poudin , to ^ obtain, d by finite-difference solution 

oi tae rto.-dy case of tne air injection problem. file initial 
pr f.'tuu: u die tr ibut lOU p\ x , C ,) = p^ u ^(xj is calculated by tbe 
Caine or ueoduri ,-,o employed for the calculation of f^ and 

( V/ i 

•k but k b vii.il be replaced by k b £ . Pig. 7.1 snows the 
various bouxidary and space coxxditions for rhe problem analysed 
h.,r. i; 'll iriati cully . idle solution may be obtained by the 
nuim rio.ui methods described in Chapters 111, IV and V for the 
sturdy state problems after me time-dependent term is 
repi'tOod by taio backward finite difference approximation. 


iln, eiutier-bmith method described in Chapters III and IV 
w “s tri d first ,.nd it- failed to converge for the reasons similai 
to tuot.u given in Chapter V. l’her fore, trie finite-difference 
m thod of Chapter V was used with tne details provided in the 
next : i; ron . 


7 ,;j i u.ito Uiffu.ence Method for the Unsteady homogeneous 
j nj c c 1 1 on £ r o oi t-m 


iu ux- 


of.M.jy opuations are solved at discrete time- 


atepf! t i , i = 1,2,3, . . . with step-sices progressively 
increasing in the positiv r e t-dir action % 


t l = % 


t b - t ip-1 = At ± = k t At i _ 1 , i = 2,5, . 


(7.8a) 

(7.8b) 



on ADPB : f+4x 


on ADCB 


VPo^, x 

X7ZT7 ~^“ ( ~ 


f«=0, H=H h or H'=0 


on GCBH. : f’=l, H=1 


f=f^°\ H=H p=p^°' ) 


on AEHE : f=f » H=H 1 , p=p. 


Pig# 7.1. Diagramatic arrangement of the boundary 
and initial conditions 
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ILuiv liiu iirst time step At ^ is taken of 0(10”*^) and 
' K t; ^ ^ * 'i-’he ;irr,viigoment of t and x stations is shown below. 


,x i, Ax 2 

. '«!**» > * *4 

. , . . i 




4x j 


At 


! X 1 x ;; *7 X j-2 x 3-l x 7 *J-2 X J-1 X J =1 


i i 


->X 


5’ 



# 

» 


t. 


i-2 


At. , 

JL- J 


Ax j = k x Ax j-l» 3 ~ 

Ati = k(. 3 = 2,3,4,... 

k x , k t > 1 


*1-1 


At i 



i 

t 


Figure 7.2 


Schematic Arrangement of x and t Stations 
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(i 


' w * ,U1 *^ 1 unction <*>(x,t) the first derivative in the 
•u"i u ( I form can be written as 


, cj oi . ( -U 




u. _ ( U;L + Y ^) + v i c/ 1 " 2 ^ (7.9a) 


wh x 1 . 


/ 0 td V ( . ! ) 


l. «<3) _ (o + q ) uj^- 1 ) 


+ q . co 
D 


,(3-2) 


i “ 

m . 1 
« 'l Ax 

tJ 


V. = u 

X 


u 


for i = 1 
3 = 1,2 


.. _i_ ^V 1 ) 

i “ At i ^ Efc 

l, = -i- ("IT—) 

.1 AX . ■ K x 
3 


r 4 = -1- p— ) 

i At i '“Jfc + 1 

.2 

1 (■ X 


) 


for i > 2 

3 > 3 


(7.9b) 


(7.10a) 


(7.10b) 


3 AXj k x + l 

with the* help of finite-difference equations (5.44) and 
(7.9) equations (7.1) and (7.2), give (with F = f ' ) : 

f x f 


x f 1 


+ P r—^ + 




2(Aij) 2 ~ 4 Un) (An)J f U2 (Ati) 2 4 (An) (An) 




+ 


— ^ - 4x »j + p(y I - 1) n] 

(An) 2 b J p 
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** 1 * J) -p 2C f 

”7 ~P (i' G - 2F_ + Fj (p - p ) _ £ (P - P ) 

2 (An)" c b a 4 (An) V c a (An) c a 


ix dp 


+ k (~Z~ - 1) H - 4x F_ S_ ( jt- + q . ) + 4x J. q. • 

P C b bl 3 3 b bl 3 

> ( x- 1 ) , , T -mi— 2 


-t- 4x 


| ~(u. + y . ) F'-* + y. p 

L 1 1 e l e 


] 


(7.11) 


P 1 


x f. 


'’ L ' (An) 2 4(4 n) (Ati)-J ' "f L M 2( An) 2 ' 4(An) 




1 


+ H„ 


p i . £ 

T 


t ^ + H n f— ' ^“o - 4x *** *, - 4x u. + 4p x 


( An) -J 0 L p r ( An) 2 e 3 


i-1 i-2 

u i V - ( u i + v i^ P + v i P 


— (H - 2H_ + H ) - -4-r ( H ~ - Hj 
ii>r(An) 2 G b 


a ' 4 ( An) e a 


•.X S',, H_ ( * , + (1, ) + 4x J H~" c q. - T-r-S (H. - H = ) 


ri-2 


X f 


o 


‘j TH f e o ^3 (An) e a 


2p 

+ Fr <* - Pr) 


S' 


2( An) 


(P„ - 2»_ + S a + - 2S 0 + I? d ) 


2 v c 


b 


4 ( An) 2 
w i-2 

+ V • 

i e 


<* f - V o 0 - * a ) 


+ 4x 


'K + v i> H o _1 


te— i— 1 i-2*1 g 

4(3 S u. p - (u i + v j _) p + v x P _P e 

n L. 


(7.12) 
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.in equations (7*11) and (7.12) the variables without the 
jjupcfiiot’ij t ra fur to their values on the ith time-step and 

Mk 

vanohA with subscript ’bl’stand for their values at the 

ht -iturfu Subscripts * b * and ‘e* refer to the (j-l)th 
( j;ti* wtwti oiis respectively along the plate. If we put = v^ = 0 
ir. eis* (7*11 } uhd (7.12), the steady case equations are retrieved. 

Since wo are using backward difference formulas for the 

ti-..e-<ioi , ivr-.tivea, they result in an implicit scheme which is 

inherently it table. At each time step t = t. , 1 = 2,3,4.*,, we 

( 1 ) 

n; proximate p by extrapolation of values at the previous 
M:rj-nt.;pri w explained later. Then the system of algebraic 
..■luntiona resulting Trow (7.11) and (7.12) is solved in the 
region x, < x < x,, 0 < n < n using the initial condition 
(7,3 - b) and the boundary conditions (7.6 - 7). From the 
solution wo then obtain ^ in accordance with the interaction 
<: iuati ona (Si. 5 0) and (2.51) with s 0. The solution then 

procut-da^to the next time-step 1: = t. + 1, It is continued 
until j < 1C" 4 and is then compared with the steady state 

solution corronding to 0 * 6^ £• 

-(i) 

Tno extrapolation formula required for p is as fallows : 

p(xj - (l+ktJpfe^ - (^ 13) 

Ji) ' ' -(0) 

For the first time-step (i»l) p is taken to 

(±) . . _ 
the value for t - 0. How dp /Sx oao be obtained from p or 

can bo extrapolated from the previous time-steps using a formula 
similar to (7.13). 
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'Clio time-dependent term in eq . (7.16) be finlte-diffeneneea 


_(i) 


Ft 



( u i +v i> 


-(i-* 1 ) _(i-2 ) 

6 + v^6 


(7.13) 


When eqs. (7*18) and (7.17) are substituted into equation 
(7.16), v/e .jet; an ordinary non-linear first order differential— 
e j nation with the initial value of 6= 6(x 1 ,t) specified at 
•di tii.x— so tops . IT his can be integrated by the fourth-order 

’am u-Kh U;a method to obtain 6 and then p can be computed 
fro. >: 

p(x,t) = I(x f t )/6(x,t ) (7.19) 


1 V L t/ ii 

CO 

Kx,fc) = / [) H(x,n , t )-f ' ^(x,ri,t ) 3 dn/l ( 7 . 20 ) 

0 

where i is defined in secifion 3.3 of chaDter III. 

U 

2t may be mentioned here that equation (7; 19) is another 
form of osi . (2.50) with = 0; 


7.4 .Results of tho Unsteady Problem with and without Air 
injection 

Pig. 7.3 through 7.5 present the distribution of 
P (x,0,t), f”(x,0,t) and H'(x,0;t) with time for no-injection 

case i.e, a = 0.0 with Pr = i72, Y = 1.4, H-^ = 0.5, = 8, 

I, I = 20 rnd the change in angle of attack from 2° to 6°. The 

OQ 

initial and final states are also indicated in these figures. 

The figures also contain the location of the condition 
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b U c£/x' = 1 which correspond to the discontinuity quoted in 
the literature [~5, 3(7] . The discontinuity appears to have 
been avoided by suitable formulation of the variables, and in 
particularly avoiding the variable of the type t*u^/x*. 

lig . 7,6 through 7,8 show the variation of p(x,t), 
f ’ ’ (x,0,t) and H*(x,0,t) with time for the case of air 
injection (oc=0.15) with P = .72, y = 1.4, = .5, x L *= 8, 

J ' ;: oo = ^0 011 & "bhe change in angle of attach from 2° to 6°. 

The initial and final steady states are also indicated. Here 

even though the value of p(x,t) is smaller compared with the 

no-in joction case, the value of p would be higher because of 
the somewhat larger value of p Q in the relation: p. = (p 0 x) p. 
Pig, 7.9 contains the variation of 6 along the plate for different 
times for air injection (a=.l5) with Pr = .72, Y = 1.4, H^= «5, 

X-j^* = 8, 1^= 20 and the change in angle of attach from 2° to 6°, 

In Pig. 7.10 through 7-. 12 are given the distribution of. 
p(x,fc), f ' ’ (x, 0, t) and H* (x,0, t ) with time for the case with 
suction at the wall (a = -0.5) with Pr = .72, Y = 1,4, H^= 0.5, 
X^* = 8, M w = 20, and the change in angle of attach from 2° to 

6°. Tho initial and final steady states are also shown in 

these figures. The value of p calculated from p for this 
case would bo lower as compared to the no-injection case due 1 

to the smaller value of p Q in the relation: p = (p Q x) p. 
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Finally, Fig. 7.13 gives the variation of 6 along the 
plate for different times for the suction case (a = -.5) with 
Pr = .72, y = 1.4, = .5, X L * = 8, M m = 20 and the change 

in angle of attach from 2° to 6°. 

Figures 7.3 through 7.13 indicate that the induced 
pressure, the wall shea^ and heat transfer functions and the 
boundary layer thickness approach the final steady state values 
<i t t ^ 10. Closer to the leading edge, the steady state is 
reached earlier. 





Fig. 7. 5 Variation of H'(o) with time for P r =.72 , =0-0 , H b = .5 , 

= 20 y X* = 8 y change of angle from 2° to 6°(+ value att=0 
x final steady state 7 • location of t*u£,/x* = 1) 
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CHAPTER. Till 


COhCLUblhe REUflEKS 

Tile problem of high Mach number flow past a wedge has been 
solved, Tne special features of the present work are 

1. Hie solution for tne entire region x _> x^» where x^ 
represents the start of tne continuum regime, has been obtained 

2. The starting profiles at x=x^ (x^=G(10“^J) have been 
obtained by series solution and thus the error of entirely 
neglecting tne non-continuum processes in the region 

u < x < x 1 n«s been avoided, 

3. ho use of similarity is made and hence it was not necossar^ 
to find solutions separately for tne strong, transition and 
weak interaction sub-regions. 

4. The effect of surface mass transfix — both injection and 
suction — has 'to von consider ed by suitably modifying the 
boundary conditions at the surface. As no use of similarity 
is mr'.dc, the rate of surface mass transfer need not be 
subjected to any particular law regarding x-dependence. The 
study mciuJes the homogeneous case (the gas injected is the 
same as that in tno atmosphere) and the heterogeneous case 
(wnen tne injected gas is different from that in the 
atmosphere) . 

5. The study also includes considerations of the problem of 

the. time t aten for the boundary layer to adjust from the 
state corresponding to tne plate inclination 0 ^ t° that 

cor re so ending to 6, P when the inclination of the plate is 
impulsively cEumged from the value ®h ? i 1:10 y b,f* 
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-Che m\in results obtr.in.ed may be summed up as follows. 

1* tn o surface mass transfer reduces tne heat transfer 
appreciably. '-Che results are shown in fig. (5.23 ) 

2. Che liit and drag coefficients, the lift to drag ratio, 
the position of tne aerodynamic centre and the heat transfer 
rate per unit widtn of the wedge-wing have been calculated 
and they -ro presented in the form of a table on page I5u. 

3. iW pressure distribution for various values of a — the 
injec tion/suction parameter — has been snown in fig. (3.8; 

4 . i'nt change in the value of shin friction coefficient 
uni boun lary layer thickness due to injection is shown 
respectively in figs. (5. 22) and (5.24), This data is 

pr^" ents 1 lor throe different gases - helium, air and argon, 
b. i‘or tne study of the problem of the adjustment of the 
boundary layer when the inclination of the plate is impulsively 
chang. '1 from B b • to G b ^ -fche variation from the initial 
to find values of pressure, sh^ar stress function, wall 
jnthalpy gradient and boundary layer thickness is shown for 
thro*- c •"> s ; s : u, = U, a > U (injection; and a < 0 (suction) 
in fi, ur* s (7.3) to (7.13; . 

lutur* work u:,y be directed towards 
(i) obfc.-ining tne series solutions for a binary gas mixture. 
( n ) solving the time-dependent problem with heterogeneous 
injection 

(iiij case (li; with gradual change in the .angle of attack. 
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APPENDIX A 


DEl'AlJjfc) OP SIMPLIFICATIONS USED IE 
U lif Aih I k G 1’iIE FLOW GOVEKHING EQUATIONS 


lxi this appendix we will describe the details of 
simplifications of some of the terms required in the non- 
d imonsionalisation of the equations (2.2) and (2.4)* The 
various terms which require these details are 


i is! Lf aja J 1>A' 

7 * ta* ’ e* at* 1=1 1 ay* 

i'irwt we simplify •% • il: ° m th ® e< i uation o£ state 

p 3 x 


or 


«* = R* P* T* 


* R * 

m, Poo _* _ u 

S' *4 p 1 with TT* " R a " 


• R. 


T 

H oo oo 


m 


A 


the Universal gas constant and R* is the gas 
constant for the mixture. 
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whore 0 * = f ( 0 . 0 * ) lias boon used. 

P £ J- P X 

Furthur, for a binary mixture i = 1,2 and 

C 1 + c 2 = 1 

Hero o x is for the injected gas and c g is for the main-stream 
gas, namely, air. Therefore we car write i 
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If wo substitute for p from (A.l) and eliminate I from the 
resulting expression with the help of (A. 2) we can write, 
nit or minor simplifications, 
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APPENDIX B 


MODIFIED SIMPSON'S RULE 

Pile integrals S^ (i = 1,2, 3, 4) of Chapter YI require a 
special lormula for evaluation if the step-size along the plate 
i,j not uniform. In the case of heterogeneous injection, where 
tno results are obtained with the finite-difference method, 
tnore is no difficulty in the evaluation of these integrals 
since the x-stations employed in the computations are at 
equal intervals and hence Simpson's rule can be conveniently 
uscjd. But in the case of Clutter-Smith method a non-uniform 
stop-size given by ax. = k ax. 1 is used to reduce the 
computational time. Therefore, a modification is to be made 
in the Simpson's rule of integration to use it for the non- 
uniform space -intervals. 

Let 0 1 , Og, Og,..., 0^ ... Oj represent the value of 


the integrands of any of the integrals S^ (i = 1,2,3,4). Then 
the evaluation of the integrals can be made as follows : 
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i+e lolling to Fig. B«l, let the 0x00, between the ordinates 
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one middle point of x^ and x^ +1 can be obtained as 


° 1+ 1 - 0. 
0 . = 0 • + — «i= — . .. — J . 

3, mean J k AT 

2C AX 3 


AX. (1 + k ) 


AX 


i] (Ba) 


•rfith u. >UWQn given by (B.l), we can get the value of A^ using 


the Simpson' f.'. rule*: 
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i ; ovj , putting 3 = 2 ,4, 6 , 8 , .. . (J-l) and summing up all the 
resulting values of A. we get, the value of the 
integral (i = 1,2, 3, 4 ) as : 
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*For the Simpson's rule with equal intervals, A^ is given by : 
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v/hore ax^ = Ax g has been used and J is an odd integer. 

By putting ^ = 1 in (B,3) we can obtain Simpson’s rule 
Tor equal intervals, Cn in eq. (B.3) corresponds to the value 
of fcho respective integrands appearing in the integrals S^^ 

(i = 1, 2,3,4). 

The values of the integrals S. (i = 1, 2,3,4) obtained 
by using (B.3) was checked in three different ways % 

(i) by fitting a cubic spline through the upper ends 
of the ordinates and applying Simpson's rule to obtain the 
area after computing the ordinates at equal intervals from the 
cubic spline. 

(ii) by the trapezoidal rule applied on each segment 
bounded by two adjacent ordinates. 

and (iii) by plotting a graph and then measuring the area. 

In the first of the above three methods of checking, 
the results are found to agree upto fourth places after the 
decimal point; arid with the other two methods a good agreement 
of the results has been noted. 

As mentioned earlier, J is an odd integer as it appears 
in oq. (B.3), Bor an even value of J, expression (B.3) can be 
used upto (J-l)th station and the area of the last segment may 
be obtained by the trapezoidal rule. 



ilPPMIX C 


UTILJ^TIOU OP SERIES EXPANSION SOLUTIONS 
IV. I Til AND tflTUOUT AIR INJECTION 


'L'ho expressions contained in this section for the 
evaluation of various aerodynamic coefficients are similar 
to those obtained in ref. £7l3» In 'the evaluation of these 
coefficients for the case with air injection, however, 
appropriate solutions to the boundary layer equations obtained 
in 0 nap ter IJi must be used. 


Tin; wall shear stress t* can bo expressed in the form 
3/2 b 

T b 2 * 


where the constants'^ (i = 0, 1,2,3) are given by 
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Substituting for p* and t* from the expressions (3.8) and (G.l) 
tne equations (6.7) and (6.8) yield, after necessary integrations 
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whore 


- = i \*AL 


Prom (0.3) and (0.4) we may also write 


where 


and 
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ihe subscripts p and v denote the contributions of the pressure 
and viscous forces. 


Pig, 0.1 compares tile presently obtained 1/1 ratio for 
a flat plate (e^ = 0) with the results of Mir els and 1 ewe Hen 
£ 72[] and Rodkiewicz and Ohattopadhyay Q71]] . A good agree- 
ment is seen for small values of angle of attack. Por large 
values of angle of attack, present computations are closer to 



C3 


the result;;; of Miroic and Lcwollen £72]] . It may be pointed 
,,ut, hur>; that tiho difference between the present results and 
those obtained. by Rodkiewicz and Chattopadhyay arises due to 
uie variations in the values of the various constants obtained 
in the two calculations. As mentioned in Sec. 3.6 of Chapter 
ill, Rodkiewicz and Chattopadhyay have evaluated their 
constants by integrating the governing equations near the wall 
with a Rungc-Kutta method employing a large step-size (An = 0.10). 
jn our calculations, however, a more precise method, namely, 
tho Taylor’s series expansion with much smaller step-size 
(0.10/16) was used to initiate the solutions at the wall. 







